Chapter Two 



Some Concepts 
Before Bioreactors 



INTRODUCTION 



B 



efore understanding bioreactors we must understand some basic 
concepts of fluid dynamics the properties and parameters of fluids 
and the mediums used in the bioreactors. 



Fluid Dynamics 



Fluid dynamics is a sub-discipline of fluid mechanics that deals with fluid 
flow— the natural science of fluids (liquids and gases) in motion. It has several 
sub-disciplines itself, including aerodynamics (the study of air and other gases 
in motion) and hydrodynamics (the study of liquids in motion). Fluid dynamics 
has a wide range of applications, but here we will only discuss its application on 
bioreactors. Some of its principles are even used in traffic engineering, where 
traffic is treated as a continuous fluid. 

Fluid dynamics offers a syste-matic structure — which underlies these 
practical disciplines— that embraces empirical and semi-empirical laws derived 
from flow measurement and used to solve practical problems. The solution to a 
fluid dynamics problem typically involves calculating various properties of the 
fluid, such as velocity, pressure, density, and temperature, as functions of space 
and time. 
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Historically, hydrodynamics meant something different than it does today. 
Before the twentieth century, hydrodynamics was synonymous with fluid 
dynamics. This is still reflected in names of some fluid dynamics topics, like 
magneto-hydrodynamics and hydrodynamic stability — both also applicable in, 
as well as being applied to, gases. 

Equations of Fluid Dynamics 

The foundational axioms of fluid dynamics are the conservation laws, 
specifically, conservation of mass, conservation of linear momentum (also 
known as Newton's Second Law of Motion), and conservation of energy 
(also known as First Law of Thermodynamics). These are based on classical 
mechanics and are modified in quantum mechanics and general relativity. They 
are expressed using the Reynolds Transport Theorem. In addition to the above, 
fluids are assumed to obey the continuum assumption. Fluids are composed 
of molecules that collide with one another and solid objects. However, the 
continuum assumption considers fluids to be continuous, rather than discrete. 
Consequently, properties such as density, pressure, temperature, and velocity 
are taken to be well-defined at infinitesimally small points, and are assumed to 
vary continuously from one point to another. The fact that the fluid is made up 
of discrete molecules is ignored. For fluids which are sufficiently dense to be a 
continuum, do not contain ionized species, and have velocities small in relation 
to the speed of light, the momentum equations for Newtonian fluids are the 
Navier-Stokes equations, which is a non-linear set of differential equations that 
describes the flow of a fluid whose stress depends linearly on velocity gradients 
and pressure. The unsimplified equations do not have a general closed-form 
solution, so they are primarily of use in Computational Fluid Dynamics. The 
equations can be simplified in a number of ways, all of which make them easier 
to solve. Some of them allow appropriate fluid dynamics problems to be solved 
in closed form. In addition to the mass, momentum, and energy conservation 
equations, a thermodynamical equation of state giving the pressure as a function 
of other thermodynamic variables for the fluid is required to completely specify 
the problem. An example of this would be the perfect gas equation of state: 

P M 

where p is pressure, p is density, R u is the gas constant, M is the molar mass and 
T is temperature. 
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Compressible vs Incompressible Flow 



All fluids are compressible to some extent, that is, changes in pressure or 
temperature will result in changes in density. However, in many situations the 
changes in pressure and temperature are sufficiently small that the changes in 
density are negligible. In this case the flow can be modeled as an incompressible 
flow. Otherwise the more general compressible flow equations must be used. 

Mathematically, incompressibility is expressed by saying that the density p 
of a fluid parcel does not change as it moves in the flow field, i.e., 

D* 

where D/Dt is the substantial derivative, which is the sum of local and convective 
derivatives. This additional constraint simplifies the governing equations, 
especially in the case when the fluid has a uniform density. For flow of gases, to 
determine whether to use compressible or incompressible fluid dynamics, the 
Mach number of the flow is to be evaluated. As a rough guide, compressible 
effects can be ignored at Mach numbers below approximately 0.3. For liquids, 
whether the incompressible assumption is valid depends on the fluid properties 
(specifically the critical pressure and temperature of the fluid) and the flow 
conditions (how close to the critical pressure the actual flow pressure becomes). 
Acoustic problems always require allowing compressibility, since sound waves 
are compression waves involving changes in pressure and density of the medium 
through which they propagate. 

Viscous vs Inviscid Flow 

Viscous problems are those in which fluid friction has significant effects on 
the fluid motion. 

The Reynolds number, which is a ratio between inertial and viscous forces, 
can be used to evaluate whether viscous or inviscid equations are appropriate 
to the problem. 

Stokes flow is flow at very low Reynolds numbers, Re<<l, such that inertial 
forces can be neglected compared to viscous forces. On the contrary, high 
Reynolds numbers indicate that the inertial forces are more significant than the 
viscous (friction) forces. Therefore, we may assume the flow to be an inviscid 
flow, an approximation in which we neglect viscosity completely, compared to 




Fig. 3: Potential flow around a wing 
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inertial terms. This idea can work fairly 

well when the Reynolds number is 

high. However, certain problems such 

as those involving solid boundaries, . 

may require that the viscosity be 

included. Viscosity often cannot 

be neglected near solid boundaries 

because the no-slip condition can 

generate a thin region of large strain 

rate (known as Boundary layer) 

which enhances the effect of even a 

small amount of viscosity, and thus 

generating vorticity. Therefore, to calculate net forces on bodies (such as wings) 

we should use viscous flow equations. As illustrated by d'Alembert s paradox, a 

body in an inviscid fluid will experience no drag force. The standard equations 

of inviscid flow are the Euler equations. Another often used model, especially in 

computational fluid dynamics, is to use the Euler equations away from the body 

and the boundary layer equations, which incorporates viscosity, in a region 

close to the body. The Euler equations can be integrated along a streamline to 

get Bernoulli's equation. When the flow is everywhere irrotational and inviscid, 

Bernoulli's equation can be used throughout the flow field. Such flows are called 

potential flows. 

Steady vs Unsteady Flow 

When all the time derivatives of a flow field vanish, the flow is considered 
to be a steady flow. Steady-state flow refers to the condition where the fluid 
properties at a point in the system do not change over time. Otherwise, flow is 
called unsteady. Whether a particular flow is steady or unsteady, can depend 
on the chosen frame of reference. For instance, laminar flow over a sphere is 
steady in the frame of reference that is stationary with respect to the sphere. 
In a frame of reference that is stationary with respect to a background flow, the 
flow is unsteady. Turbulent flows are unsteady by definition. A turbulent flow 
can, however, be statistically stationary. According to Pope: The random field 
U(x,t) is statistically stationary if all statistics are invariant under a shift in time. 
This roughly means that all statistical properties are constant in time. Often, 
the mean field is the object of interest, and this is constant too in a statistically 
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Fig. 4: Hydrodynamics simulation of the Rayleigh-Taylor instability 
stationary flow. Steady flows are often more tractable than otherwise similar 
unsteady flows. The governing equations of a steady problem have one dimension 
fewer (time) than the governing equations of the same problem without taking 
advantage of the steadiness of the flow field. 

Laminar vs Turbulent Flow 

Turbulence is flow characterized by recirculation, eddies, and apparent 
randomness. Flow in which turbulence is not exhibited is called laminar. It 
should be noted, however, that the presence of eddies or recirculation alone 
does not necessarily indicate turbulent flow — these phenomena may be present 
in laminar flow as well. Mathematically, turbulent flow is often represented via 
a Reynolds decomposition, in which the flow is broken down into the sum of an 
average component and a perturbation component. It is believed that turbulent 
flows can be described well through the use of the Navier-Stokes equations. 
Direct numerical simulation (DNS), based on the Navier-Stokes equations, 
makes it possible to simulate turbulent flows at moderate Reynolds numbers. 
Restrictions depend on the power of the computer used and the efficiency of 
the solution algorithm. The results of DNS have been found to agree well with 
experimental data for some flows. Most flows of interest have Reynolds numbers 
much too high for DNS to be a viable option, given the state of computational 
power for the next few decades. Any flight vehicle large enough to carry a 
human (L > 3 m), moving faster than 72 km/h (20 m/s) is well beyond the limit 
of DNS simulation (Re = 4 million). Transport aircraft wings (such as on an 
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Airbus A300 or Boeing 747) have Reynolds numbers of 40 million (based on the 
wing chord). In order to solve these real-life flow problems, turbulence models 
will be a necessity for the foreseeable future. Reynolds-averaged Navier-Stokes 
equations (RANS) combined with turbulence modeling provides a model of the 
effects of the turbulent flow. Such a modeling mainly provides the additional 
momentum transfer by the Reynolds stresses, although the turbulence also 
enhances the heat and mass transfer. Another promising methodology is large 
eddy simulation (LES), especially in the guise of detached eddy simulation 
(DES) — which is a combination of RANS turbulence modeling and large eddy 
simulation. 

Newtonian vs Non-Newtonian Fluids 



Sir Isaac Newton showed how stress and the rate of strain are very close to 
linearly related for many familiar fluids, such as water and air. These Newtonian 
fluids are modeled by a coefficient called viscosity, which depends on the specific 
fluid. However, some of the other materials, such as emulsions and slurries and 
some visco-elastic materials (e.g. blood, some polymers), have more complicated 
non-Newtonian stress-strain behaviours. These materials include sticky liquids 
such as latex, honey, and lubricants which are studied in the sub-discipline of 
rheology. 

Subsonic vs Transonic, Supersonic and Hypersonic Flows 

While many terrestrial flows (e.g. flow of water through a pipe) occur at low 
mach numbers, many flows of practical interest (e.g. in aerodynamics) occur at 
high fractions of the Mach Number M=l or in excess of it (supersonic flows). 
New phenomena occur at these Mach number regimes (e.g. shock waves for 
supersonic flow, transonic instability in a regime of flows with M nearly equal 
to 1, non-equilibrium chemical behavior due to ionization in hypersonic flows) 
and it is necessary to treat each of these flow regimes separately. 

Magnetohydrodynamics 

Magnetohydrodynamics is the multi-disciplinary study of the flow of 
electrically conducting fluids in electromagnetic fields. Examples of such fluids 
include plasmas, liquid metals, and salt water. The fluid flow equations are 
solved simultaneously with Maxwells equations of electromagnetism. This is 
necessary because we will see in the future during the design of the bioreactor 
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when the electrodes and sensors of different instruments are introduced in the 
bioreactor vessel, the Eddy Current is developed, which need to be carefully 
tackled by using a good knowledge of Magnetothermodynamics. 

Other Approximations 

There are a large number of other possible approximations to fluid dynamic 
problems. Some of the more commonly used are listed below. 

• The Boussine sq approximation neglects variations in density except 
to calculate buoyancy forces. It is often used in free convection 
problems where density changes are small. 

• Lubrication theory and Hele-Shaw flow exploits the large aspect ratio 
of the domain to show that certain terms in the equations are small 
and so can be neglected. 

• Slender-body theory (regarding the vessel of the bioreactor) is a 
methodology used in Stokes flow problems to estimate the force on, 
or flow field around, a long slender object in a viscous fluid( used and 
agar medium in bioreactors). 

• The shallow-water equations can be used to describe a layer of 
relatively inviscid fluid with a free surface, in which surface gradients 
are small. 

• The Boussinesq equations are applicable to surface waves on thicker 
layers of fluid and with steeper surface slopes. 

• Darcy s law is used for flow in porous media, and works with variables 
averaged over several pore-widths. 

• In rotating systems, like bioreactors the quasi-geostrophic 
approximation assumes an almost perfect balance between pressure 
gradients and the Coriolis force. It is useful in the study of atmospheric 
dynamics. 

Terminology in Fluid Dynamics 

The concept of pressure is central to the study of both fluid statics and 
fluid dynamics. A pressure can be identified for every point in a body of fluid, 
regardless of whether the fluid is in motion or not. Pressure can be measured 
using an aneroid, Bourdon tube, mercury column, or various other methods. 
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Some of the terminology that is necessary in the study of fluid dynamics is not 
found in other similar areas of study. In particular, some of the terminology 
used in fluid dynamics is not used in fluid statics. 

Terminology in Incompressible Fluid Dynamics 

The concepts of total pressure and dynamic pressure arise from Bernoulli's 
equation and are significant in the study of all fluid flows. (These two pressures 
are not pressures in the usual sense — they cannot be measured using an 
aneroid, Bourdon tube or mercury column.) To avoid potential ambiguity 
when referring to pressure in fluid dynamics, many authors use the term static 
pressure to distinguish it from total pressure and dynamic pressure. Static 
pressure is identical to pressure and can be identified for every point in a fluid 
flow field. In Aerodynamics, L.J. Clancy writes: To distinguish it from the total 
and dynamic pressures, the actual pressure of the fluid, which is associated not 
with its motion but with its state, is often referred to as the static pressure, but 
where the term pressure alone is used it refers to this static pressure. A point in 
a fluid flow where the flow has come to rest {i.e. speed is equal to zero adjacent 
to some solid body immersed in the fluid flow) is of special significance. It is of 
such importance that it is given a special name— a stagnation point. The static 
pressure at the stagnation point is of special significance and is given its own 
name — stagnation pressure. In incompressible flows, the stagnation pressure at 
a stagnation point is equal to the total pressure throughout the flow field. 

Terminology in Compressible Fluid Dynamics 

In a compressible fluid, such as air, the temperature and density are 
essential when determining the state of the fluid. In addition to the concept 
of total pressure (also known as stagnation pressure), the concepts of total (or 
stagnation) temperature and total (or stagnation) density are also essential in any 
study of compressible fluid flows. To avoid potential ambiguity when referring 
to temperature and density, many authors use the terms static temperature and 
static density. Static temperature is identical to temperature; and static density 
is identical to density; and both can be identified for every point in a fluid flow 
field. The temperature and density at a stagnation point are called stagnation 
temperature and stagnation density. A similar approach is also taken with the 
thermodynamic properties of compressible fluids. Many authors use the terms 
total (or stagnation) enthalpy and total (or stagnation) entropy. The terms static 
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enthalpy and static entropy appear to be less common, but where they are 
used they mean nothing more than enthalpy and entropy respectively, and the 
prefix "static" is being used to avoid ambiguity with their 'total' or stagnation 
counterparts. Because the 'total' flow conditions are defined by isentropically 
bringing the fluid to rest, the total (or stagnation) entropy is by definition always 
equal to the "static" entropy. 

FLUIDIZED BED 

A fluidized bed is formed when a quantity of a solid particulate substance 
(usually present in a holding vessel) is placed under appropriate conditions to 
cause the solid/fluid mixture to behave as a fluid. This is usually achieved by the 
introduction of pressurized fluid through the particulate medium. This results 
in the medium then having many properties and characteristics of normal 
fluids; such as the ability to free-flow under gravity, or to be pumped using fluid 
type technologies. 

The resulting phenomenon is called fluidization. Fluidized beds are used 
for several purposes, such as fluidized bed reactors (types of chemical reactors), 
fluid catalytic cracking, fluidized bed combustion, heat or mass transfer 
or interface modification, such as applying a coating onto solid items. This 
technique is also becoming more common in Aquaculture for the production 
of shellfish in Integrated Multi-Trophic Aquaculture systems. 

History 

In 1922 Fritz Winkler made the first industrial application of fluidization in 
a reactor for a coal gasification process. In 1942, the first circulating fluid bed was 
built for catalytic cracking of mineral oils, with fluidization technology applied 
to metallurgical processing (roasting arsenopyrite) in the late 1940s. During 
this time theoretical and experimental research improved the design of the 
fluidized bed. In the 1960s VAW-Lippewerk in Lunen, Germany implemented 
the first industrial bed for the combustion of coal and later for the calcination 
of aluminium hydroxide. 

Properties of Fluidized Beds 

A fluidized bed consists of fluid-solid mixture that exhibits fluid-like 
properties. As such, the upper surface of the bed is relatively horizontal, which 
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is analogous to hydrostatic behavior. The bed can be considered to be an 
inhomogeneous mixture of fluid and solid that can be represented by a single 
bulk density. Furthermore, an object with a higher density than the bed will sink, 
whereas an object with a lower density than the bed will float, thus the bed can 
be considered to exhibit the fluid behavior expected of Archimedes' principle. 
As the "density", (actually the solid volume fraction of the suspension), of the 
bed can be altered by changing the fluid fraction, objects with different densities 
comparative to the bed can, by altering either the fluid or solid fraction, be 
caused to sink or float. In fluidized beds, the contact of the solid particles with 
the fluidization medium (a gas or a liquid) is greatly enhanced when compared 
to packed beds. This behavior in fluidized combustion beds enables good 
thermal transport inside the system and good heat transfer between the bed 
and its container. Similarly to the good heat transfer, which enables thermal 
uniformity analogous to that of a well mixed gas, the bed can have a significant 
heat-capacity whilst maintaining a homogeneous temperature field. 
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Fig. 4: A diagram of a fluidized bed 
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Fluidized Bed Types 



Bed types can be coarsely classified by their flow behavior, including : 

• Stationary or bubbling bed is the classical approach where the gas 
at low velocities is used and fluidization of the solids is relatively 
stationary, with some fine particles being entrained. 

• Circulating Fluidized Beds (CFB), where gases are at a higher velocity 
sufficient to suspend the particle bed, due to a larger kinetic energy 
of the fluid. As such the surface of the bed is less smooth and larger 
particles can be entrained from the bed than for stationary beds. 
Entrained particles are recirculated via an external loop back into the 
reactor bed. Depending on the process, the particles may be classified 
by a cyclone separator and separated from or returned to the bed, 
based upon particle cut size. 

• Vibratory Fluidized beds are similar to stationary beds, but add 
a mechanical vibration to further excite the particles for increased 
entrainment. 

• Transport or flash reactor (FR). At velocities higher than CFB, 
particles approach the velocity of the gas. Slip velocity between gas 
and solid is significantly reduced at the cost of less homogeneous heat 
distribution. 

• Annular fluidized bed (AFB). A large nozzle at the center of a bubble 
bed introduces gas as high velocity achieving the rapid mixing zone 
above the surrounding bed comparable to that found in the external 
loop of a CFB. 

Basic Model 

When the packed bed has a fluid passed over it, the pressure drop of the 
fluid is approximately proportional to the fluid s superficial velocity. In order 
to transition from a packed bed to a fluidized condition, the gas velocity is 
continually raised. For a free-standing bed there will exist a point, known as the 
minimum or incipient fluidization point, whereby the bed s mass is suspended 
directly by the flow of the fluid stream. The corresponding fluid velocity, known 
as the "minimum fluidization velocity" Wm/. Beyond the minimum fluidization 
velocity ( u > u^f ), the bed material will be suspended by the gas-stream and 
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further increases in the velocity will have a reduced effect on the pressure, 
owing to sufficient percolation of the gas flow. Thus the pressure drop from for 
u > Umf is relatively constant. At the base of the vessel the apparent pressure 
drop multiplied by the cross-section area of the bed can be equated to the force 
of the weight of the solid particles (less the buoyancy of the solid in the fluid). 

Ap w = H W (1- e w ){p 5 - p f )g 
Geldart Groupings 



In 1973, Professor D. Geldart proposed the grouping of powders in to 
four so-called "Geldart Groups". The groups are defined by their locations on a 
diagram of solid-fluid density difference and particle size. Design methods for 
fluidized beds can be tailored based upon the particle s Geldart grouping: 

Group A: For this group the particle size is between 20 and 100 ^m, and 
the particle density is typically less than 1.4g/cm 3 . Prior to the initiation of a 
bubbling bed phase, beds from these particles will expand by a factor of 2 to 3 at 
incipient fluidization, due to a decreased bulk density. Most powder-catalyzed 
beds utilize this group. 

Group B: The particle size lies between 40 and 500 ^m and the particle 
density between 1.4-4g/cm 3 . Bubbling typically forms directly at incipient 
fluidization. 

Group C: This group contains extremely fine and consequently the most 
cohesive particles. With a size of 20 to 30 um, these particles fluidize under very 
difficult to achieve conditions, and may require the application of an external 
force, such as mechanical agitation. 

Group D: The particles in this region are above 600 ^m and typically have 
high particle densities. Fluidization of this group requires very high fluid energies 
and is typically associated with high levels of abrasion. Drying grains and peas, 
roasting coffee beans, gasifying coals, and some roasting metal ores are such 
solids, and they are usually processed in shallow beds or in the spouting mode. 

Distributor 

Typically, pressurized gas or liquid enters the fluidized bed vessel through 
numerous holes via a plate known as a distributor plate, located at the bottom 
of the fluidized bed. The fluid flows upward through the bed, causing the solid 
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particles to be suspended. If the inlet fluid is disabled the bed may settle or pack 
onto the plate. 

Application 

Fluidized beds are used as a technical process which has the ability to 
promote high levels of contact between gases and solids. In a fluidized bed a 
characteristic set of basic properties can be utilized, indispensable to modern 
process and chemical engineering, these properties include: 

• Extremely high surface area contact between fluid and solid per unit 
bed volume. 

• High relative velocities between the fluid and the dispersed solid 
phase. 

• High levels of intermixing of the particulate phase. 

• Frequent particle-particle and particle-wall collisions. 

In the Application of Bioreactors: In the bioreactors the Agar and other 
nutrients medium is also applied in the form of a bed in the bioreactor vessel 
where the bacteria and also antibiotics are cultured. 

Taking an example from the food processing industry: fluidized beds are 
used to accelerate freezing in some IQF tunnel freezers. IQF means Individually 
Quick Frozen, or freezing unpackaged separate pieces. These fluidized bed 
tunnels are typically used on small food products like peas, shrimp or sliced 
vegetables, and may use cryogenic or vapor-compression refrigeration. The 
fluid used in fluidized beds may also contain a fluid of catalytic type; that's why 
it is also used to catalyze the chemical reaction and also to improve the rate of 
reaction. 

SOME STANDARD LAWS AND 
EQUATION GOVERNING THE FLUID 
DYNAMICS IN BIOREACTOR 

Fluid management in bioreactor, some mathematical laws are applied to 
manipulate the processes taking place inside the reactor. There are now being 
explained, with their applications as well. 



Bioreaclor: lis Fundamentals, Design and Applications 



Darcy-Weisbach Equation 



In fluid dynamics, the Darcy-Weisbach equation is a phenomenological 
equation, which relates the head loss — or pressure loss — due to friction along 
a given length of pipe to the average velocity of the fluid flow. The equation is 
named after Henry Darcy and Julius Weisbach. 

The Darcy-Weisbach equation contains a dimensionless friction factor, 
known as the Darcy friction factor. This is also called the Darcy-Weisbach 
friction factor or Moody friction factor. The Darcy friction factor is four times 
the Fanning friction factor, with which it should not be confused. 

History 

Historically this equation arose as a variant on the Prony equation; this 
variant was developed by Henry Darcy of France, and further refined into the 
form used today by Julius Weisbach of Saxonyin 1845. Initially, data on the 
variation of f with velocity was lacking, so the Darcy-Weisbach equation was 
outperformed at first by the empirical Prony equation in many cases. In later 
years it was eschewed in many special-case situations in favor of a variety of 
empirical equations valid only for certain flow regimes, notably the Hazen- 
Williams equation or the Manning equation, most of which were significantly 
easier to use in calculations. However, since the advent of the calculator, ease 
of calculation is no longer a major issue, and so the Darcy-Weisbach equations 
generality has made it the preferred one. 

Derivation 

The Darcy-Weisbach equation is a phenomenological formula obtainable 
by dimensional analysis. Away from the ends of the pipe, the characteristics of 
the flow are independent of the position along the pipe. The key quantities are 
then the pressure drop along the pipe per unit length, Ap/L, and the volumetric 
flow rate. The flow rate can be converted to an average velocity V by dividing by 
the wetted area of the flow (which equals the cross-sectional area of the pipe if 
the pipe is full of fluid). 

Pressure has dimensions of energy per unit volume. Therefore, the 
pressure drop between two points must be proportional to (l/2)pV 2 y which has 
the same dimensions as it resembles (see below) the expression for the kinetic 
energy per unit volume. We also know that pressure must be proportional to 
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the length of the pipe between the two points L as the pressure drop per unit 
length is a constant. To turn the relationship into a proportionality coefficient 
of dimensionless quantity we can divide by the hydraulic diameter of the pipe, 
D, which is also constant along the pipe. Therefore, 

Ap oc £ - \ P V\ 

The proportionality coefficient is the dimensionless "Darcy friction factor" 
or "flow coefficient". This dimensionless coefficient will be a combination of 
geometric factors such as tt, the Reynolds number and (outside the laminar 
regime) the relative roughness of the pipe (the ratio of the roughness height to 
the hydraulic diameter). 

Note that (1/2) pV 2 is not the kinetic energy of the fluid per unit volume, 
for the following reasons. Even in the case of laminar flow, where all 
the flow lines are parallel to the length of the pipe, the velocity of the 
fluid on the inner surface of the pipe is zero due to viscosity, and the 
velocity in the center of the pipe must therefore be larger than the 
average velocity obtained by dividing the volumetric flow rate by the 
wet area. The average kinetic energy then involves the mean-square 
velocity, which always exceeds the square of the mean velocity. In the 
case of turbulent flow, the fluid acquires random velocity components 
in all directions, including perpendicular to the length of the pipe, and 
thus turbulence contributes to the kinetic energy per unit volume but 
not to the average lengthwise velocity of the fluid. 

Head Loss Form 

Head loss can be calculated with 

where 

hAs the head loss due to friction (SI units: m); 
L is the length of the pipe (m); 

D is the hydraulic diameter of the pipe (for a pipe of circular section, this 
equals the internal diameter of the pipe) (m); 
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V is the average velocity of the fluid flow, equal to the volumetric flow 

rate per unit cross-sectional wetted area (m/s); 
gis the local acceleration due to gravity (m/s 2 ); 
f D is a dimensionless coefficient called the Darcy friction factor. It can 

be found from a Moody diagram or more precisely by solving the 

Colebrook equation. Do not confuse this with the Fanning Friction 

factor,/ 

Pressure Loss Form 

Given that the head loss by expresses the pressure loss Ap as the height of 
a column of fluid, 

Ap = p ■ g ■ k f 

where p is the density of the fluid, the Darcy- Weisbach equation can also be 
written in terms of pressure loss: 

where the pressure loss due to friction Ap (units: Pa or kg/ms 2 ) is a function 
of: 

the ratio of the length to diameter of the pipe, L/D; 

the density of the fluid, p (kg/m 3 ); 

the mean velocity of the flow, V (m/s), as defined above; 

Darcy Friction Factor; a (dimensionless) coefficient of laminar, or 

turbulent flow, f D . 
Since the pressure loss equation can be derived from the head loss 
equation by multiplying each side by p and g. 

Darcy Friction Factor 

The friction factor f D or flow coefficient A is not a constant and depends 
on the parameters of the pipe and the velocity of the fluid flow, but it is known 
to high accuracy within certain flow regimes. It may be evaluated for given 
conditions by the use of various empirical or theoretical relations, or it may be 
obtained from published charts. These charts are often referred to as Moody 
diagrams, after L. F. Moody, and hence the factor itself is sometimes called the 
Moody friction factor. It is also sometimes called the Blasius friction factor, 
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after the approximate formula he proposed. For laminar (slow) flows, it is a 
consequence of Poiseuille s law that \ = 64/Re, where Re is the Reynolds number 
calculated substituting for the characteristic length the hydraulic diameter of 
the pipe, which equals the inside diameter for circular pipe geometries. For 
turbulent flow, methods for finding the friction factor f include using a diagram 
such as the Moody chart; or solving equations such as the Colebrook- White 
equation, or the Swamee-Jain equation. While the diagram and Colebrook- 
White equation are iterative methods, the Swamee-Jain equation allows/to be 
found directly for full flow in a circular pipe. 

Confusion with the Fanning Friction Factor 

The Darcy-Weisbach friction factor, f D is 4 times larger than the Fanning 
friction factor, f, so attention must be paid to note which one of these is meant 
in any "friction factor" chart or equation being used. Of the two, the Darcy- 
Weisbach factor, f D is more commonly used by civil and mechanical engineers, 
and the Fanning factor, f, by chemical engineers, but care should be taken to 
identify the correct factor regardless of the source of the chart or formula. 

Note that 

L pV 2 L a 

Most charts or tables indicate the type of friction factor, or at least provide 
the formula for the friction factor with laminar flow. If the formula for laminar 
flow is f = 16/Re, its the Fanning factor, f, and if the formula for laminar flow is 
f D = 64/Re, its the Darcy-Weisbach factor, f D . 

Which friction factor is plotted in a Moody diagram may be determined by 
inspection if the publisher did not include the formula described above: 

Observe the value of the friction factor for laminar flow at a Reynolds 
number of 1000. 

If the value of the friction factor is 0.064, then the Darcy friction factor is 
plotted in the Moody diagram. Note that the nonzero digits in 0.064 are the 
numerator in the formula for the laminar Darcy friction factor: f D = 64/Re. 

If the value of the friction factor is 0.016, then the Fanning friction factor 
is plotted in the Moody diagram. Note that the nonzero digits in 0.016 are the 
numerator in the formula for the laminar Fanning friction factor: f = 16/Re. 
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The procedure above is similar for any available Reynolds number that is 
an integral power often. It is not necessary to remember the value 1000 for this 
procedure— only that an integral power of ten is of interest for this purpose. 

Darcy Friction Factor Formulae 

The Darcy friction factor formulae are equations — based on experimental 
data and theory — for the Darcy friction factor. The Darcy friction factor 
is a dimensionless quantity used in the Darcy- Weisbach equation, for the 
description of friction losses in pipe flow as well as open channel flow. It is also 
known as the Darcy- Weisbach friction factor or Moody friction factor and is 
four times larger than the Fanning friction factor. 

Flow Regime 

Which friction factor formula may be applicable depends upon the type of 
flow that exists: 

• Laminar Flow 

• Transition between laminar and turbulent flow 

• Fully turbulent flow in smooth conduits 

• Fully turbulent flow in rough conduits 

• Free surface flow. 

LAMINAR FLOW 

The Darcy friction factor for laminar flow (Reynolds number less than 
2000) is given by the following formula: 

3 Re 
where: 

/ is the Darcy friction factor 

Re is the Reynolds number. 

TRANSITION FLOW 

Transition (neither fully laminar nor fully turbulent) flow occurs in the 
range of Reynolds numbers between 2300 and 4000. The value of the Darcy 
friction factor may be subject to large uncertainties in this flow regime. 
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TURBULENT FLOW IN SMOOTH CONDUITS 

Empirical correlations exist for this flow regime. Such correlations are 
included in the ASHRAE Handbook of Fundamentals. 

TURBULENT FLOW IN ROUGH CONDUITS 

The Darcy friction factor for fully turbulent flow (Reynolds number greater 
than 4000) in rough conduits is given by the Colebrook equation. 

FREE SURFACE FLOW 

The last formula in the Colebrook equation section of this article is for free 
surface flow. The approximations elsewhere in this article are not applicable for 
this type of flow. 

Choosing a Formula 

Before choosing a formula it is worth knowing that in the paper on the 
Moody chart, Moody stated the accuracy is about ±5% for smooth pipes and 
±10% for rough pipes. If more than one formula is applicable in the flow regime 
under consideration, the choice of formula may be influenced by one or more 
of the following: 

• Required precision 

• Speed of computation required 

• Available computational technology: 

—calculator (minimize keystrokes) 
— spreadsheet (single-cell formula) 
— programming/scripting language (subroutine). 
Colebrook Equation (Compact Forms) 

The Colebrook equation is an implicit equation that combines experimental 
results of studies of turbulent flow in smooth and rough pipes. It was developed 
in 1939 by C. F. Colebrook. The 1937 paper by C. F. Colebrook and C. M. White 
is often erroneously cited as the source of the equation. This is partly because 
Colebrook in a footnote (from his 1939 paper) acknowledges his debt to White 
for suggesting the mathematical method by which the smooth and rough 
pipe correlations could be combined. The equation is used to iteratively solve 
for the Darcy- Weisbach friction factor f. This equation is also known as the 
Colebrook- White equation. 
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For conduits that are flowing completely full of fluid at Reynolds numbers 
greater than 4000, it is defined as: 



1 / c 2.51 \ 

7j = - 2h ^{37D h + ^77) 



or 

1 _. / e 2.51 \ 

where: 

/ is the Darcy friction factor 

Roughness height, 6 (m, ft) 

Hydraulic diameter, D^ (m, ft) — For fluid-filled, circular conduits, Dh = 

D = inside diameter 
Hydraulic radius, 7? h (m, ft) — For fluid-filled, circular conduits, 7? h = D/4 

= (inside diameter) /4 
R e is the Reynolds number. 
Solving 

The Colebrook equation used to be solved numerically due to its apparent 
implicit nature. Recently, the Lambert W function has been employed to 
obtained explicit reformulation of the Colebrook equation. 
Expanded Forms 

Additional, mathematically equivalent forms of the Colebrook equation 
are: 



77 = ir384 - 2,os '»U + W7) 



where: 

1.7384... = 2 log (2 x 3.7) = 2 log (7.4) 
18.574 = 2.51 x 3.7 x 2 
and 



-L = 1 .l 3 64... + 21o gloW£ )-21o glo (l + S j^_) 
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or 



where: 

1.1364... = 1.7384... - 2 log (2) = 2 log (7.4) - 2 log (2) = 2 log (3.7) 
9.287 = 18.574/2 = 2.51 x 3.7. 

The additional equivalent forms above assume that the constants 3.7 
and 2.51 in the formula at the top of this section are exact. The constants are 
probably values which were rounded by Colebrook during his curve fitting; but 
they are effectively treated as exact when comparing (to several decimal places) 
results from explicit formulae (such as those found elsewhere in this article) to 
the friction factor computed via Colebrook s implicit equation. 

Equations similar to the additional forms above (with the constants 
rounded to fewer decimal places — or perhaps shifted slightly to minimize 
overall rounding errors) may be found in various references. It may be helpful 
to note that they are essentially the same equation. 

Free Surface Flow 

Another form of the Colebrook- White equation exists for free surfaces. 
Such a condition may exist in a pipe that is flowing partially full of fluid. For free 
surface flow: 

1 / e 2.51 \ 

77 = " 21ogl0 li2^ + ^77J- 

Approximations of the Colebrook equation 

Haaland Equation 

The Haaland equation is used to solve directly for the Darcy-Weisbach 
friction factor / for a full-flowing circular pipe. It is an approximation of the 
implicit Colebrook- White equation, but the discrepancy from experimental 
data is well within the accuracy of the data. It was developed by S. E. Haaland 
in 1983. 

The Haaland equation is defined as: 
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where: 

/is the Darcy friction factor 

e/Dis the relative roughness 

Re is the Reynolds number. 
Swamee-Jain Equation 

The Swamee-Jain equation is used to solve directly for the Darcy- Weisbach 
friction factor / for a full-flowing circular pipe. It is an approximation of the 
implicit Colebrook- White equation. 

0.25 



/ = 



Niod^+S^)]' 



where/is a function of: 

Roughness height, 8 (m, ft) 

Pipe diameter, D (m, ft) 

Reynolds number, Re (unitless). 
Serghides's Solution 

Serghides's solution is used to solve directly for the Darcy- Weisbach 
friction factor / for a full-flowing circular pipe. It is an approximation of the 
implicit Colebrook- White equation. It was derived using Steffensens method. 

The solution involves calculating three intermediate values and then 
substituting those values into a final equation. 



,4 = -2 log 



B = -2 log. 



10 



(e/D 12 \ 

UjD 2.5L4 \ 
^3.7 + Re j 



.'e/D 2.515 
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f= y A ~ C-2B + A) 



where/is a function of: 

Roughness height, e (m, ft) 

Pipe diameter, D (m, ft) 

Reynolds number, Re (unitless). 

The equation was found to match the Colebro ok- White equation within 
0.0023% for a test set with a 70-point matrix consisting often relative roughness 
values (in the range 0.00004 to 0.05) by seven Reynolds numbers (2500 to 108). 

Goudar-Sonnad Equation 

Goudar equation is the most accurate approximation to solve directly for 
the Darcy-Weisbach friction factor f for a full-flowing circular pipe. It is an 
approximation of the implicit Colebrook- White equation. Equation has the 
following form: 

2 
*"b(10) 

EfD 

ln(10)ife 
d = ^02~ 

s = bd + In(d) 



q = s s/(s+1) 

d 
g = bd + In - 

z = In — 



D r .A = Z 
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Dcfa = D LA ( 



ff+1 

1 + 



*/2 



( 3 +l) 2 + (z/3)(2 5 -l), 
1 



7f = a 



["I D t , 



where/is a function of: 

Roughness height, 8 (m, ft) 

Pipe diameter, D (m, ft) 

Reynolds number, Re (unitless). 
Brki Solution 

Brkic shows one approximation of the Colebrook equation based on the 
Lambert W- function: 

S = ln Re 



l-8161ns^a^ 



1 



where Darcy friction factor f is a function of: 

Roughness height, 8 (m, ft) 

Pipe diameter, D (m, ft) 

Reynolds number, Re (unitless). 

The equation was found to match the Colebrook- White equation within 
3.15%. 
Blasius Correlations 

Early approximations by Blasius are given in terms of the Fanning friction 
factor in the Paul Richard Heinrich Blasius article. 

Table of Approximations 

The following table lists historical approximations where: 

Re, Reynolds number (unitless) 
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A, Darcy friction factor (dimensionless) 

£, roughness of the inner surface of the pipe (m, ft) 

D, inner pipe diameter (m, ft) 

KOZENY-CARMAN EQUATION 

The equation was first proposed by Kozeny (1927) and later modified by 
Carman (1937, 1956). 

The Kozeny-Carman equation (or Carman-Kozeny equation) is a relation 
used in the field of fluid dynamics to calculate the pressure drop of a fluid 
flowing through a packed bed of solids. It is named after Josef Kozeny and Philip 
C. Carman. The equation is only valid for laminar flow. 

The equation is given as: 



Ap 18QV^(l-£) 



2 



where Ap is the pressure drop, L is the total height of the bed, Vq I s the superficial 
or "empty- tower" velocity, (A is the viscosity of the fluid, e is the porosity of the bed, 
O s is the sphericity of the particles in the packed bed, andD is the diameter of the 
related spherical particle. This equation holds for flow through packed beds with 
particle Reynolds numbers up to approximately 1.0, after which point frequent 
shifting of flow channels in the bed causes considerable kinetic energy losses. 
This equation can be expressed as "flow is proportional to the pressure drop and 
inversely proportional to the fluid viscosity" which is known as Darcy s law. 

Reynolds Number 

Reynolds number (Re) is a dimensionless number that gives a measure 
of the ratio of inertial forces (which characterize how much a particular fluid 
resists any change in motion, not to be confused with inertial forces defined 
in the classical way) to viscous forces and consequently quantifies the relative 
importance of these two types offerees for given flow conditions. 

The concept was introduced by George Gabriel Stokes in 1851, but 
the Reynolds number is named after Osborne Reynolds (1842-1912), who 
popularized its use in 1883. 
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Reynolds numbers frequently arise when performing dimensional analysis 
of fluid dynamics problems, and as such can be used to determine dynamic 
similitude between different experimental cases. 

They are also used to characterize different flow regimes, such as laminar 
or turbulent flow: laminar flow occurs at low Reynolds numbers, where 
viscous forces are dominant, and is characterized by smooth, constant fluid 
motion; turbulent flow occurs at high Reynolds numbers and is dominated by 
inertial forces, which tend to produce chaotic eddies, vortices and other flow 
instabilities. 

History 



Osborne Reynolds FRS (23 August 1842 - 21 February 1912) was a 
prominent innovator in the understanding of fluid dynamics. 

Definition 

Reynolds number can be defined for a number of different situations where 
a fluid is in relative motion to a surface. These definitions generally include the 
fluid properties of density and viscosity, plus a velocity and a characteristic 
length or characteristic dimension. This dimension is a matter of convention — 
for example a radius or diameter are equally valid for spheres or circles, but 
one is chosen by convention. For aircraft or ships, the length or width can be 
used. For flow in a pipe or a sphere moving in a fluid the internal diameter is 
generally used today. Other shapes such as rectangular pipes or non-spherical 
objects have an equivalent diameter defined. For fluids of variable density such 
as compressible gases or fluids of variable viscosity such non-Newtonian fluids, 
special rules apply. The velocity may also be a matter of convention in some 
circumstances, notably stirred vessels. 

pvL vL 

Ke = = — 

/.* v 

where: 

\^is the mean velocity of the object relative to the fluid (SI units: m/s) 
I is a characteristic linear dimension, (travelled length of the fluid; 

hydraulic diameter when dealing with river systems) (m) 
^ is the dynamic viscosity of the fluid (Pa.s or N.s/m 2 or kg/(m.s)) 
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v is the kinematic viscosity (V = f-t/p) (m 2 /s) 
p is the density of the fluid (kg/m 3 ) 

Note that multiplying the Reynolds number by j- yields — g-, which is the 
ratio of the inertial forces to the viscous forces. It could also be considered the 
ratio of the total momentum transfer to the molecular momentum transfer. 

Flow in Pipe 



For flow in a pipe or tube, the Reynolds number is generally defined as: 

Re = P vD ^ = yD h = QDh 
// v vA 

where: 

D H is the hydraulic diameter of the pipe; its characteristic travelled 

length, L, (m). 
Q is the volumetric flow rate (m 3 /s). 
A is the pipe cross-sectional area (m 2 ). 

v is the mean velocity of the object relative to the fluid (SI units: m/s). 
(A is the dynamic viscosity of the fluid (Pa*s or N*s/m 2 or kg/(m*s)). 
v is the kinematic viscosity (v = t^/p) (m 2 /s). 
p is the density of the fluid (kg/m 3 ). 

For shapes such as squares, rectangular or annular ducts where the height 
and width are comparable, the characteristic dimension for internal flow 
situations is taken to be the hydraulic diameter, D H , defined as: 

44 

B p 

where A is the cross-sectional area and P is the wetted perimeter. The wetted 
perimeter for a channel is the total perimeter of all channel walls that are in 
contact with the flow. This means the length of the water exposed to air is not 
included in the wetted perimeter. 

For a circular pipe, the hydraulic diameter is exactly equal to the inside 
pipe diameter, as can be shown mathematically. 

For an annular duct, such as the outer channel in a tube-in-tube heat 
exchanger, the hydraulic diameter can be shown algebraically to reduce to 
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where 

D Q is the inside diameter of the outside pipe, and 

D i is the outside diameter of the inside pipe. 

For calculations involving flow in non-circular ducts, the hydraulic 
diameter can be substituted for the diameter of a circular duct, with reasonable 
accuracy. 

Flow in a Wide Duct 



For a fluid moving between two plane parallel surfaces— where the width 
is much greater than the space between the plates — then the characteristic 
dimension is twice the distance between the plates. 

Flow in an Open Channel 

For flow of liquid with a free surface, the hydraulic radius must be determined. 
This is the cross-sectional area of the channel divided by the wetted perimeter. 
For a semi-circular channel, it is half the radius. For a rectangular channel, the 
hydraulic radius is the cross-sectional area divided by the wetted perimeter. 
Some texts then use a characteristic dimension that is four times the hydraulic 
radius, chosen because it gives the same value of Re for the onset of turbulence 
as in pipe flow, while others use the hydraulic radius as the characteristic length- 
scale with consequently different values of Re for transition and turbulent flow. 

Object in a Fluid 

The Reynolds number for an object in a fluid, called the particle Reynolds 
number and often denoted Re , is important when considering the nature of 
flow around that grain, whether or not vortex shedding will occur, and its fall 
velocity. 

Sphere in a Fluid 

Note that purely laminar flow 



For a sphere in a fluid, the 
characteristic length-scale is the diameter 
of the sphere and the characteristic 
velocity is that of the sphere relative 
to the fluid some distance away from 
the sphere, such that the motion of the 



only exists up to Re = 0. 1 under this 
definition. Under the condition of 
low Re, the relationship between 
force and speed of motion is given 
by Stokes' law. 



Some Concepts Before Bioreactor 

sphere does not disturb that reference parcel of fluid. The density and viscosity 
are those belonging to the fluid. 

Oblong Object in a Fluid 



The equation for an oblong object is identical to that of a sphere, with the 
object being approximated as an ellipsoid and the axis of length being chosen 
as the characteristic length scale. Such considerations are important in natural 
streams, for example, where there are few perfectly spherical grains. For grains 
in which measurement of each axis is impractical, sieve diameters are used 
instead as the characteristic particle length-scale. Both approximations alter the 
values of the critical Reynolds number. 

Fall Velocity 



The particle Reynolds number is important in determining the fall velocity 
of a particle. When the particle Reynolds number indicates laminar flow, Stokes' 
law can be used to calculate its fall velocity. When the particle Reynolds number 
indicates turbulent flow, a turbulent drag law must be constructed to model the 
appropriate settling velocity. 

Packed Bed 



For fluid flow through a bed of approximately spherical particles of diameter 
D in contact, if the voidage is e and the superficial velocity is V, the Reynolds 
number can be defined as: 



Laminar conditions apply up to Re = 10, fully turbulent from 2000. 
Stirred Vessel 



In a cylindrical vessel stirred by a central rotating paddle, turbine or 
propeller, the characteristic dimension is the diameter of the agitator D. The 
velocity is ND where N is the rotational speed. Then the Reynolds number is: 

pND 2 
Re = - . 

The system is fully turbulent for values of Re above 10000. 
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Transition and Turbulent Flow 



In boundary layer flow over a flat plate, experiments confirm that, after 
a certain length of flow, a laminar boundary layer will become unstable and 
turbulent. This instability occurs across different scales and with different fluids, 
usually when Re^, ft* 5 X 1Q 5 > where x is the distance from the leading edge of 
the flat plate, and the flow velocity is the freestream velocity of the fluid outside 
the boundary layer. 

For flow in a pipe of diameter D, experimental observations show that for 
"fully developed" flow, laminar flow occurs when Re D < 2300 and turbulent 
flow occurs when Re £? > 4000. In the interval between 2300 and 4000, laminar 
and turbulent flows are possible and are called "transition" flows, depending 
on other factors, such as pipe roughness and flow uniformity. This result is 
generalized to non-circular channels using the hydraulic diameter, allowing a 
transition Reynolds number to be calculated for other shapes of channel. These 
transition Reynolds numbers are also called critical Reynolds numbers, and 
were studied by Osborne Reynolds around 1895. 

Pipe Friction 

Pressure drops seen for fully developed flow of fluids through pipes can be 
predicted using the Moody diagram which plots the Darcy-Weisbach friction 
factor /against Reynolds number Re and relative roughnesa e/D. The diagram 
clearly shows the laminar, transition, and turbulent flow regimes as Reynolds 

M-.Mrik ] }i;i".r;Liit 
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Fig 5: The Moody diagram, which describes the Darcy-Weisbach friction factor/as a function of 
the Reynolds number and relative pipe roughness. 
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number increases. The nature of pipe flow is strongly dependent on whether the 
flow is laminar or turbulent. 

Similarity of Flows 



In order for two flows to be similar they must have the same geometry, 
and have equal Reynolds numbers and Euler numbers. When comparing fluid 
behavior at corresponding points in a model and a full-scale flow, the following 
holds: 

Re m = Re 

P™ P 



Eu m = Eu i.e. 



p™v m 2 pv 2 ' 



quantities marked with 'm' concern the flow around the model and the others 
the actual flow. This allows engineers to perform experiments with reduced 
models in water channels or wind tunnels, and correlate the data to the actual 
flows, saving on costs during experimentation and on lab time. Note that true 
dynamic similitude may require matching other dimensionless numbers as well, 
such as the Mach number used in compressible flows, or the Froude number 
that governs open-channel flows. Some flows involve more dimensionless 
parameters than can be practically satisfied with the available apparatus and 
fluids, so one is forced to decide which parameters are most important. For 
experimental flow modeling to be useful, it requires a fair amount of experience 
and judgement of the engineer. 

Smallest Scales of Turbulent Motion 

In a turbulent flow, there is a range of scales of the time-varying fluid motion. 
The size of the largest scales of fluid motion (sometimes called eddies) are set 
by the overall geometry of the flow. For instance, in an industrial smoke stack, 
the largest scales of fluid motion are as big as the diameter of the stack itself. 
The size of the smallest scales is set by the Reynolds number. As the Reynolds 
number increases, smaller and smaller scales of the flow are visible. In a smoke 
stack, the smoke may appear to have many very small velocity perturbations or 
eddies, in addition to large bulky eddies. In this sense, the Reynolds number is 
an indicator of the range of scales in the flow. The higher the Reynolds number, 
the greater the range of scales. The largest eddies will always be the same size; 
the smallest eddies are determined by the Reynolds number. 
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What is the explanation for this phenomenon? 
A large Reynolds number indicates that viscous 
forces are not important at large scales of the flow. 
With a strong predominance of inertial forces over 
viscous forces, the largest scales of fluid motion 
are undamped— there is not enough viscosity to 
dissipate their motions. The kinetic energy must 
"cascade" from these large scales to progressively 
smaller scales until a level is reached for which 
the scale is small enough for viscosity to become 
important (that is, viscous forces become of the 
order of inertial ones). It is at these small scales 
where the dissipation of energy by viscous action 
finally takes place. The Reynolds number indicates at what scale this viscous 
dissipation occurs. Therefore, since the largest eddies are dictated by the flow 
geometry and the smallest scales are dictated by the viscosity, the Reynolds 
number can be understood as the ratio of the largest scales of the turbulent 
motion to the smallest scales. 

In physiology 




Fig. 6: Creeping flow past a 

sphere: streamlines, drag force 

F d and force by gravity F . 



Poiseuille s law on blood circulation in the body is dependent on laminar 
flow. In turbulent flow the flow rate is proportional to the square root of the 
pressure gradient, as opposed to its direct proportionality to pressure gradient 
in laminar flow. Using the definition of the Reynolds number we can see that 
a large diameter with rapid flow, where the density of the blood is high, tends 
towards turbulence. Rapid changes in vessel diameter may lead to turbulent 
flow, for instance when a narrower vessel widens to a larger one. Furthermore, a 
bulge of atheroma may be the cause of turbulent flow, where audible turbulence 
may be detected with a stethoscope. 

In Viscous Fluids 

Where the viscosity is naturally high, such as polymer solutions and 
polymer melts, flow is normally laminar. The Reynolds number is very small 
and Stokes' Law can be used to measure the viscosity of the fluid. Spheres are 
allowed to fall through the fluid and they reach the terminal velocity quickly, 
from which the viscosity can be determined. 
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The laminar flow of polymer solutions is exploited by animals such as fish 
and dolphins, who exude viscous solutions from their skin to aid flow over 
their bodies while swimming. It has been used in yacht racing by owners who 
want to gain a speed advantage by pumping a polymer solution such as low 
molecular weight polyoxyethylene in water, over the wetted surface of the hull 
It is, however, a problem for mixing of polymers, because turbulence is needed 
to distribute fine filler (for example) through the material. Inventions such as 
the "cavity transfer mixer" have been developed to produce multiple folds into 
a moving melt so as to improve mixing efficiency. The device can be fitted onto 
extruders to aid mixing. 

Derivation 

The Reynolds number can be obtained when one uses the nondimensional 
form of the incompressible Navier-Stokes equations: 



(£ + H 



-VjD + juV 2 v + f 



Each term in the above equation has the units of a "body force" (force per 
unit volume) or, equivalently, an acceleration times a density. Each term is thus 
dependent on the exact measurements of a flow. When one renders the equation 
nondimensional, that is when we multiply it by a factor with inverse units of the 
base equation, we obtain a form which does not depend directly on the physical 
sizes. One possible way to obtain a nondimensional equation is to multiply the 
whole equation by the following factor: 

D 



pV* 

where: 

\^is the mean velocity, v or v, relative to the fluid (m/s). 
D is the characteristic length, Z, (m). 
p is the fluid density (kg/m 3 ) 
If we now set: 

we can rewrite the Navier-Stokes equation without dimensions: 
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g _ 1 

where the term : D y - Rg 

Finally, dropping the primes for ease of reading: 
f\ 1 

at Re 

This is why mathematically all flows with the same Reynolds number are 
comparable. Notice also, in the above equation, as: Re — v oo the viscous terms 
vanish. Thus, high Reynolds number flows are approximately inviscid in the 
free-stream. 

Reynolds' Transport Theorem 



Also known as the Leibniz-Reynolds' transport theorem, or in short 
Reynolds theorem, is a three-dimensional generalization of the Leibniz integral 
rule which is also known as differentiation under the integral sign. The theorem 
is named after Osborne Reynolds (1842-1912). It is used to recast derivatives 
of integrated quantities and is useful in formulating the basic equations of 
continuum mechanics. 

Consider integrating f = f (x, t) over the time-dependent region fi(t) that 
has boundary dQ{t), then taking the derivative with respect to time: 

d 



d*/ii 



f d V . 

o(t) 

If we wish to move the derivative under the integral sign there are two 
issues: the time dependence of/ and the introduction of and removal of space 
from n due to its dynamic boundary. Reynolds' transport theorem provides the 
necessary framework. 

General Form 



Reynolds' transport theorem, derived in, is: 

ttJ fdV= / § dV+ / (**-=) fdA 
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in which n(x, t) is the outward-pointing unit-normal, x is a point in the 
region and is the variable of integration, dV and dA are volume and surface 
elements at X, and v b (x, t) is the velocity of the area element - so not necessarily 
the flow velocity. The function / may be vector or scalar valued. Note that 
the integral on the left hand side is a function solely of time, and so the total 
derivative has been used. 

Form for a Material Element 



In continuum mechanics this theorem is often used for material elements, 
which are parcels of fluids or solids which no material enters or leaves. If fi(f ) 
is a material element then there is a velocity function v = v(x, t) and the 
boundary elements obey 



v 6 ■ n = v ■ n. 



This condition may be substituted to obtain 



( f f dv) = / § dV + / (v ■ n)f dA . 
\Jn(t) / Jn(t) v*> Jdnit) 



Proof for a Material Element 



Let f^o be reference configuration of the region il[t). Let the motion and the 
deformation gradient be given by 

x= V CM; =* F{XJ) = V <p. 

Let */(X, t) = det[F[X. 7 t)]. Then, integrals in the current and the 
reference configurations are related by 

/ f(x,*) dV= / f[ v CM,t] A^J) dV = f f{XJ) J{XJ) dv . 

That this derivation is for a material element is implicit in the time constancy 
of the reference configuration: it is constant in material coordinates. The time 
derivative of an integral over a volume is defined as 

* (L f(x ' () dv ) = &s GL, f( *<< + a,) dv - L f(x '° dv ) ■ 

Converting into integrals over the reference configuration, we get 

|- ( f ffx, ") dv) = lim -J- f / ffX,f + At) J(X,t ■+■ At) dV c - / f (X, t) J(X,t) dV^ . 
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Since STq is independent of time, we have 



& Vftffl / Ato 



f(X, (f -h At) J{X 7 t + At)-f{X,t) J[X,t) 



ir-rfJ At 



<iv fl 



= / ^[r(x : /)./(x,f)j f ]v lh 

Now, the time derivative of rfet F is given by 

^^ = |<<fctJ?) = (detFJCV-v) = J(X,t) V ^(X,*),*) = ^PM) V-vf^t) . 
Therefore, 
'KjLi/** ' )llV ) = /* (<?I [ ' (X '' )1 ■ , ( X -0 + f(X.0./(X.i) V-v(x,l)J ,1V, 
= J (|(f(X,*)| + f (X, t) V ■ v(x,*)) J(X,«) dV 

- / ffM + f^Ov-vfx,^ dv 

J(HO v ' 

where f is the material time derivative off. Now, the material derivative is 
given by 

r(x,t) = ^^+[Vf(x,t)]-v(x,t). 
Therefore, 

£(L r(x,i)tlv )~L (^t^ + i^^oj-v^o+^o v-v(x,^ ,iv 

or, 



^(L fdV )=L(^ +Vfv+fVv ) dV 



Using the identity 

V ■ (v ® w) = v(V ■ w) + Vv ■ w 

we then have 



^(i/ dv H,>(l +v( H dv - 

Using the divergence theorem and the identity (a ® b) ■ n = (b ■ n) we have 
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~(f fdv)=/ %M+[ (f®v)-nclA=/ ?^V+/ [vn)fdA d 

Erroneous Sources 



This theorem is widely quoted, incorrectly, as being the form specific 
to material volumes. See the planetmath external link below for an example. 
Clearly, if the material volume form is applied to regions other than material 
volumes, errors will ensue. 

A Special Sase 



If we take £1 to be constant with respect to time, then V^ = and the 
identity reduces to 



<HJb J a dt 



as expected. This simplification is not possible if an incorrect form of the 
Reynolds transport theorem is used. 

Interpretation and Reduction to One Dimension 

The theorem is the higher dimensional extension of Differentiation under 
the integral sign and should reduce to that expression in some cases. Suppose/ 
is independent of y & z, and that S7(t) is a unit square in the y-z plane and has 
x limits a(t) and &(£). Then Reynolds transport theorem reduces to 

which is the expression given on Differentiation under the integral sign, except 
that there the variables x and t have been swapped. 

REYNOLDS DECOMPOSITION 

Reynolds decomposition is a mathematical technique to separate the 
average and fluctuating parts of a quantity. For example, for a quantity u the 
decomposition would be 



u(x,y,z,t) = u{x,y,z) + u{x,y,z,t) 
where u denotes the time average of u (often called the steady component), and 
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U the fluctuating part (or perturbations). The perturbations are defined such 
that their time average equals zero. 

This allows us to simplify the Navier-Stokes equations by substituting in 
the sum of the steady component and perturbations to the velocity profile and 
taking the mean value. The resulting equation contains a nonlinear term known 
as the Reynolds stresses which gives rise to turbulence. 

Reynold's Stress 

Reynolds stress (or, the Reynolds stress tensor) is the stress tensor in a fluid 
obtained from the averaging operation over the Navier-Stokes equations to 
account for turbulent fluctuations in fluid momentum. 

Definition 



For a homogeneous fluid and an incompressible flow, the flow velocities are 
split into a mean part and a fluctuating part using Reynolds decomposition: 

Ui = ^+*4 

with u(x, t) being the flow velocity vector having components u i in the x^ 
coordinate direction (with x i denoting the components of the coordinate 
vector X). The mean velocities are determined by either time averaging, spatial 
averaging or ensemble averaging, depending on the flow under study. Further 
u- denotes the fluctuating (turbulence) part of the velocity. 

The components t\- of the Reynolds stress tensor are defined as: 



with p the fluid density, taken to be non-fluctuating for this homogeneous 
fluid. 

Another - often used - definition, for constant density, of the Reynolds 
stress components is: 



which has the dimensions of velocity squared, instead of stress. 
Averaging and the Reynolds Stress 



To illustrate, Cartesian vector index notation is used. For simplicity, 
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consider an incompressible fluid: 

Given the fluid velocity w . as a function of position and time, write the average 
fluid velocity as Up and the velocity fluctuation is u { . Then i£, z = 1J^ + ti ■ . 

The conventional ensemble rules of averaging are that 
a = a, 



a + h = a + 6, 

a-6 = aft. 

One splits the Euler equations or the Navier-Stokes equations into an 
average and a fluctuating part. One finds that upon averaging the fluid equations, 
a stress on the right hand side appears of the form pri^ufj . This is the Reynolds 
stress, conventionally written Riji 



The divergence of this stress is the force density on the fluid due to the 
turbulent fluctuations. 

Reynolds Averaging of the Navier-Stokes Equations 



For instance, for an incompressible, viscous, Newtonian fluid, the continuity 
and momentum equations — the incompressible Navier-Stokes equations — can 
be written as 

dxi 
and 

Dui dp / d 2 Ui \ 

Dt dxi \dxjdxj) 

where D/Dt is the Lagrangian derivative or the substantial derivative, 

D 3 d 

= ^7 + % * 



Dt~ dt^ J dxj- 
Defining the flow variables above with a time-averaged component and a 
fluctuating component, the continuity and momentum equations become 

fl(W + ttj) =Q 
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and 



at 



+ fe+«5) 



d{ui+ «■) 



dxj 



djp + p') 
dxi 



+ P 



d 2 (u^ + u' i ) 



dxjdxj 



Examining one of the terms on the left hand side of the momentum 
equation, it is seen that 

1% + U *J dx . ~ dx _ Wi+Vi) dx _ , 

where the last term on the right hand side vanishes as a result of the continuity 
equation. Accordingly, the momentum equation becomes 



djuj + u't) d{u i +i4)(u j +^ j ) 



Of 



dx j 



9(P + P r ) 



+ /* 



d 2 (ui+<4) 



dx-jdxj 



Now the continuity and momentum equations will be averaged. The 
ensemble rules of averaging need to be employed, keeping in mind that the 
average of products of fluctuating quantities will not in general vanish. After 
averaging, the continuity and momentum equations become 



&ui 
dx t 



= 0, 



and 



dui duiuj 

1 i 

dt dxj dxj 






dp d 7 v~ 



dx 



dxjdxj' 



Using the chain rule on one of the terms of the left hand side, it is revealed 
that 

duluj dui duj 

dxj J dxj t dxj ' 

where the last term on the right hand side vanishes as a result of the averaged 
continuity equation. The averaged momentum equation now becomes, after a 
rearrangement: 



dul dui dp d 

dt 3 dxj\ dxi dxj 



du z 



: — -P^j) i 



dxj 
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where the Reynolds stresses, pu£uj, are collected with the viscous normal and 



shear stress terms, ._ <**** 
Discussion 



The question then is, what is the value of the Reynolds stress? This has 
been the subject of intense modeling and interest, for roughly the past century. 
The problem is recognized as a closure problem, akin to the problem of closure 
in the BBGKY hierarchy. A transport equation for the Reynolds stress may be 
found by taking the outer product of the fluid equations for the fluctuating 
velocity, with itself. 

One finds that the transport equation for the Reynolds stress includes 
terms w ith higher-order correlations (specifically, the triple correlation 
l^lfrl^) as well as correlations with pressure fluctuations (i.e. momentum 
carried by sound waves). A common solution is to model these terms by simple 
ad-hoc prescriptions. It should also be noted that the theory of the Reynolds 
stress is quite analogous to the kinetic theory of gases, and indeed the stress 
tensor in a fluid at a point may be seen to be the ensemble average of the stress 
due to the thermal velocities of molecules at a given point in a fluid. Thus, by 
analogy, the Reynolds stress is sometimes thought of as consisting of an isotropic 
pressure part, termed the turbulent pressure, and an off-diagonal part which 
may be thought of as an effective turbulent viscosity. In fact, while much effort 
has been expended in developing good models for the Reynolds stress in a fluid, 
as a practical matter, when solving the fluid equations using computational 
fluid dynamics, often the simplest turbulence models prove the most effective. 
One class of models, closely related to the concept of turbulent viscosity, is 
the so-called K - e model(s), based upon coupled transport equations for the 
turbulent energy density K (similar to the turbulent pressure, i.e. the trace of 
the Reynolds stress) and the turbulent dissipation rate e. Typically, the average 
is formally defined as an ensemble average as in statistical ensemble theory. 
However, as a practical matter, the average may also be thought of as a spatial 
average over some length scale, or a temporal average. Note that, while formally 
the connection between such averages is justified in equilibrium statistical 
mechanics by the ergodic theorem, the statistical mechanics of hydrodynamic 
turbulence is currently far from understood. In fact, the Reynolds stress at any 
given point in a turbulent fluid is somewhat subject to interpretation, depending 
upon how one defines the average. 
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Reynolds-averaged Navier-Stokes Equations 



The Reynolds-averaged Navier-Stokes equations (or RANS equations) are 
time-averaged equations of motion for fluid flow. The idea behind the equations 
is Reynolds decomposition, whereby an instantaneous quantity is decomposed 
into its time-averaged and fluctuating quantities, an idea first proposed by 
Osborne Reynolds. The RANS equations are primarily used to describe turbulent 
flows. These equations can be used with approximations based on knowledge of 
the properties of flow turbulence to give approximate time-averaged solutions 
to the Navier-Stokes equations. For a stationary, incompressible Newtonian 
fluid, these equations can be written in Einstein notation as: 



P U J^— =Pfi + 



dxj dxj 



\dxj 



8 



u i 



-i% + ,i(— f^l -p^u' 3 



The left hand side of this equation represents the change in mean 
momentum of fluid element owing to the unsteadiness in the mean flow and 
the convection by the mean flow. This change is balanced by the mean body 
force, the isotropic stress owing to the mean pressure field, the viscous stresses, 
and apparent stress (-puju'j owing to the fluctuating velocity field, generally 
referred to as the Reynolds stress. This nonlinear Reynolds stress term requires 
additional modeling to close the RANS equation for solving, and has led to the 
creation of many different turbulence models. The time-average operator t is a 
Reynolds operator. 

Derivation of RANS Equations 

The basic tool required for the derivation of the RANS equations from 
the instantaneous Navier-Stokes equations is the Reynolds decomposition. 
Reynolds decomposition refers to separation of the flow variable (like velocity) 
into the mean (time-averaged) component ( U) and the fluctuating component ( 
U ). Because the mean operator is a Reynolds operator, it has a set of properties. 
One of these properties is that the mean of the fluctuating quantity being equal 
to zero( u f = 0)- Thus, u(x, t) = u(x) +u(x, t), where x = (x,y,z) is the 
position vector. Some authors[2] prefer using U instead of U for the mean term 
(since an overbar is sometimes used to represent a vector). In this case, the 
fluctuating term U is represented instead by U. This is possible because the two 
terms do not appear simultaneously in the same equation. To avoid confusion, 
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the notation u U, and U will be used to represent the instantaneous, mean, 
and fluctuating terms, respectively. 

The properties of Reynolds operators are useful in the derivation of the 
RANS equations. Using these properties, the Navier-Stokes equations of 
motion, expressed in tensor notation, are (for an incompressible Newtonian 
fluid): 

ox t 

dui &u±__ f iQp flV 

dt 3 dxj 2 p dxi dxjdxj 
where fi is a vector representing external forces. 

Next, each instantaneous quantity can be split into time-averaged and 
fluctuating components, and the resulting equation time-averaged, to yield: 

§i = ° 

OXi 






dt dxj 3 dxj pdxi dxjd; 

The momentum equation can also be written as, 



J 



dui _ dui j 1 dp d 2 Ui du^u^ 

dt 3 dxj * pdxi dxjdxj dxj 

On further manipulations this yields, 

dUi _ du t - 9 r _ _ - — -T-] 
9 ~dt + pUj d^ =pf * + lhT \~P 5 *i + 2 ^ ~ P u i u 5\ 

where, £& = - \ =— -I- ^ I is the mean rate of strain tensor. 

Finally, since integration in time removes the time dependence of the 
resultant terms, the time derivative must be eliminated, leaving: 

^dotT = P ^ + d~> i p5ij + 2p ^ ij ~ pu ^n ■ 



' j 
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Reynolds' Dilatancy 



Reynolds' dilatancy is the observed tendency of a compacted granular 
material to dilate (expand in volume) as it is sheared. This occurs because the 
grains in a compacted state are interlocking and therefore do not have the 
freedom to move around one another. When stressed, a lever motion occurs 
between neighboring grains, which produces a bulk expansion of the material. 
On the other hand, when a granular material starts in a very loose state it 
may initially compact instead of dilating under shear. Reynolds' dilatancy is a 
common feature of the soils and sands studied by geotechnical engineers, and is 
a part of the broader topic of soil mechanics. It was first described scientifically 
by Osborne Reynolds (1842-1912) in 1885 and 1886. 

Ergun Equation 

The Ergun equation, derived by the Turkish chemical engineer Sabri Ergun 
in 1952, expresses the friction factor in a packed column as a function of the 
Reynolds number: 

150 

where f p and Gr p are defined as 

where: Ap is the pressure drop across the bed, 

L is the length of the bed (not the column), 

D is the equivalent spherical diameter of the packing, 

p is the density of fluid, 

\i is the dynamic viscosity of the fluid, 

V s is the superficial velocity (i.e. the velocity that the fluid would have 

through the empty tube at the same volumetric flow rate), and 
e is the void fraction of the bed (Bed porosity at any time). 

BIOFILM 

A biofilm is an aggregate of microorganisms in which cells adhere to each 
other on a surface. These adherent cells are frequently embedded within a 
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self-produced matrix of extracellular polymeric substance (EPS). Biofilm EPS, 
which is also referred to as slime (although not everything described as slime is 
a biofilm), is a polymeric conglomeration generally composed of extracellular 
DNA, proteins, and polysaccharides. Biofilms may form on living or non-living 
surfaces and can be prevalent in natural, industrial and hospital settings. The 
microbial cells growing in a biofilm are physiologically distinct from planktonic 
cells of the same organism, which, by contrast, are single-cells that may float or 
swim in a liquid medium. Microbes form a biofilm in response to many factors, 
which may include cellular recognition of specific or non-specific attachment 
sites on a surface, nutritional cues, or in some cases, by exposure of planktonic 
cells to sub-inhibitory concentrations of antibiotics. When a cell switches to the 
biofilm mode of growth, it undergoes a phenotypic shift in behavior in which 
large suites of genes are differentially regulated. 

Formation 

Formation of a biofilm begins with the attachment of free-floating 
microorganisms to a surface. These first colonists adhere to the surface initially 
through weak, reversible adhesion via van der Waals forces. If the colonists are 
not immediately separated from the surface, they can anchor themselves more 
permanently using cell adhesion structures such as pili. Some species are not 
able to attach to a surface on their own but are often able to anchor themselves 
to the matrix or directly to earlier colonists. It is during this colonization that 
the cells are able to communicate via quorum sensing using such products as 
AHL. Once colonization has begun, the biofilm grows through a combination 
of cell division and recruitment. The final stage of biofilm formation is known 
as development, and is the stage in which the biofilm is established and may 
only change in shape and size. The development of a biofilm may allow for an 
aggregate cell colony (or colonies) to be increasingly antibiotic resistant. 
Development 

There are five stages of biofilm development (see fig. 7): 

1. Initial attachment: 

2. Irreversible attachment: 

3. Maturation I: 

4. Maturation II: 

5. Dispersion: 
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Fig. 7: Five stages of biofilm development: (1) Initial attachment, (2) Irreversible attachment, (3) 
Maturation I, (4) Maturation II, and (5) Dispersion. Each stage of development in the diagram is 
paired with a photomicrograph of a developing P. aeruginosa biofilm. All photomicrographs are 

shown to same scale. 

Dispersal 



Dispersal of cells from the biofilm colony is an essential stage of the biofilm 
life cycle. Dispersal enables biofilms to spread and colonize new surfaces. 
Enzymes that degrade the biofilm extracellular matrix, such as dispersin B 
and deoxyribonuclease, may play a role in biofilm dispersal. Biofilm matrix 
degrading enzymes may be useful as anti-biofilm agents. Recent evidence has 
shown that a fatty acid messenger, cis-2-decenoic acid, is capable of inducing 
dispersion and inhibiting growth of biofilm colonies. Secreted by Pseudomonas 
aeruginosa, this compound induces cyclo heteromorphic cells in several species 
of bacteria and the yeast Candida albicans. Nitric oxide has also been shown 
to trigger the dispersal of biofilms of several bacteria species at sub-toxic 
concentrations. Nitric oxide has the potential for the treatment of patients that 
suffer from chronic infections caused by biofilms. 
Properties 

Biofilms are usually found on solid substrates submerged in or exposed to 
an aqueous solution, although they can form as floating mats on liquid surfaces 
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and also on the surface of leaves, particularly in high humidity climates. Given 
sufficient resources for growth, a biofilm will quickly grow to be macroscopic 
(visible to the naked eye). Biofilms can contain many different types of 
microorganism, e.g. bacteria, archaea, protozoa, fungi and algae; each group 
performs specialized metabolic functions. However, some organisms will form 
single-species films under certain conditions. 

Extracellular Matrix 

The biofilm is held together and protected by a matrix of excreted 
polymeric compounds called EPS. EPS is an abbreviation for either extracellular 
polymeric substance or exopolysaccharide. This matrix protects the cells within 
it and facilitates communication among them through biochemical signals. 
Some biofilms have been found to contain water channels that help distribute 
nutrients and signalling molecules. This matrix is strong enough that under 
certain conditions, biofilms can become fossilized. Bacteria living in a biofilm 
usually have significantly different properties from free-floating bacteria of the 
same species, as the dense and protected environment of the film allows them 
to cooperate and interact in various ways. One benefit of this environment is 
increased resistance to detergents and antibiotics, as the dense extracellular 
matrix and the outer layer of cells protect the interior of the community. In 
some cases antibiotic resistance can be increased a thousand fold Lateral gene 
transfer is greatly facilitated in biofilms and leads to a more stable biofilm 
structure. 

However, biofilms are not always less susceptible to antibiotics. For 
instance, the biofilm form of Pseudomonas aeruginosa has no greater resistance 
to antimicrobials than do stationary-phase planktonic cells, although when 
the biofilm is compared to logarithmic phase planktonic cells, the biofilm does 
have greater resistance to antimicrobials. This resistance to antibiotics in both 
stationary phase cells and biofilms may be due to the presence of persister cells. 

Examples 

Biofilms are ubiquitous. Nearly every species of microorganism, not only 
bacteria and archaea, have mechanisms by which they can adhere to surfaces 
and to each other. Biofilms will form on virtually every non-shedding surface in 
a non-sterile aqueous (or very humid) environment. 
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• Biofilms can be found on rocks and pebbles at the bottom of most 
streams or rivers and often form on the surface of stagnant pools of 
water. In fact, biofilms are important components of food chains in 
rivers and streams and are grazed by the aquatic invertebrates upon 
which many fish feed. 

• Biofilms can grow in the most extreme environments: from, for 
example, the extremely hot, briny waters of hot springs ranging from 
very acidic to very alkaline, to frozen glaciers. 

• In the human environment, biofilms can grow in showers very easily 
since they provide a moist and warm environment for the biofilm to 
thrive. Biofilms can form inside water and sewage pipes and cause 
clogging and corrosion. Biofilms on floors and counters can make 
sanitation difficult in food preparation areas. 

• Biofilms in cooling- or heating-water systems are known to reduce 
heat transfer. 

• Biofilms in marine engineering systems, such as pipelines of the 
offshore oil and gas industry, can lead to substantial corrosion 
problems. Corrosion is mainly due to abiotic factors; however, at least 
20% of corrosion is caused by microorganisms that are attached to the 
metal subsurface (Le., microbially- influenced corrosion). 

• Bacterial adhesion to boat hulls serves as the foundation for biofouling 
of seagoing vessels. Once a film of bacteria forms, it is easier for 
other marine organisms such as barnacles to attach. Such fouling can 
reduce maximum vessel speed by up to 20%, prolonging voyages and 
consuming fuel. Time in dry dock for refitting and repainting reduces 
the productivity of shipping assets, and the useful life of ships is 
also reduced due to corrosion and mechanical removal (scraping) of 
marine organisms from ships' hulls. 

• Biofilms can also be harnessed for constructive purposes. For example, 
many sewage treatment plants include a treatment stage in which waste 
water passes over biofilms grown on filters, which extract and digest 
organic compounds. In such biofilms, bacteria are mainly responsible 
for removal of organic matter (BOD), while protozoa and rotifers are 
mainly responsible for removal of suspended solids (SS), including 
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pathogens and other microorganisms. Slow sand filters rely on biofilm 
development in the same way to filter surface water from lake, spring 
or river sources for drinking purposes. What we regard as clean water 
is a waste material to these microcellular organisms since they are 
unable to extract any further nutrition from the purified water. 

• Biofilms can help eliminate petroleum oil from contaminated 
oceans or marine systems. The oil is eliminated by the hydrocarbon- 
degrading activities of microbial communities, in particular by a 
remarkable recently-discovered group of specialists, the so-called 
hydrocarbonoclastic bacteria (HCB). 

• Stromatolites are layered accretionary structures formed in shallow 
water by the trapping, binding and cementation of sedimentary grains 
by microbial biofilms, especially of cyanobacteria. Stromatolites 
include some of the most ancient records of life on Earth, and are still 
forming today. 

• Biofilms are present on the teeth of most animals as dental plaque, 
where they may cause tooth decay and gum disease. 

• Biofilms are found on the surface of and inside plants. They can 
either contribute to crop disease or, as in the case of nitrogen-fixing 
Rhizobium on roots, exist symbiotically with the plant. Examples 
of crop diseases related to biofilms include Citrus Canker, Pierce s 
Disease of grapes, and Bacterial Spot of plants such as peppers and 
tomatoes. 

• Biofilms are used in microbial fuel cells (MFCs) to generate electricity 
from a variety of starting materials, including complex organic waste 
and renewable biomass. 

Biofilms and Infectious Diseases 



Biofilms have been found to be involved in a wide variety of microbial 
infections in the body, by one estimate 80% of all infections. Infectious processes 
in which biofilms have been implicated include common problems such as 
urinary tract infections, catheter infections, middle-ear infections, formation 
of dental plaque, gingivitis, coating contact lenses, and less common but more 
lethal processes such as endocarditis, infections in cystic fibrosis, and infections 
of permanent indwelling devices such as joint prostheses and heart valves. 
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More recently it has been noted that bacterial biofilms may impair cutaneous 
wound healing and reduce topical antibacterial efficiency in healing or treating 
infected skin wounds. It has recently been shown that biofilms are present on 
the removed tissue of 80% of patients undergoing surgery for chronic sinusitis. 
The patients with biofilms were shown to have been denuded of cilia and goblet 
cells, unlike the controls without biofilms who had normal cilia and goblet 
cell morphology. Biofilms were also found on samples from two of 10 healthy 
controls mentioned. The species of bacteria from interoperative cultures did 
not correspond to the bacteria species in the biofilm on the respective patient s 
tissue. In other words, the cultures were negative though the bacteria were 
present. Biofilms can also be formed on the inert surfaces of implanted devices 
such as catheters, prosthetic cardiac valves and intrauterine devices. New 
staining techniques are being developed to differentiate bacterial cells growing 
in living animals, e.g. from tissues with allergy-inflammations. 

Pseudomonas aeruginosa Biofilms 

The achievements of medical care in industrialised societies are markedly 
impaired due to chronic opportunistic infections that have become increasingly 
apparent in immunocompromised patients and the aging population. Chronic 
infections remain a major challenge for the medical profession and are of 
great economic relevance because traditional antibiotic therapy is usually 
not sufficient to eradicate these infections. One major reason for persistence 
seems to be the capability of the bacteria to grow within biofilms that protects 
them from adverse environmental factors. Pseudomonas aeruginosa is not 
only an important opportunistic pathogen and causative agent of emerging 
nosocomial infections but can also be considered a model organism for the 
study of diverse bacterial mechanisms that contribute to bacterial persistence. 
In this context the elucidation of the molecular mechanisms responsible for the 
switch from planktonic growth to a biofilm phenotype and the role of inter- 
bacterial communication in persistent disease should provide new insights 
in P. aeruginosa pathogenicity, contribute to a better clinical management of 
chronically infected patients and should lead to the identification of new drug 
targets for the development of alternative anti-infective treatment strategies. 

Legionellosis 

Legionella bacteria are known to grow under certain conditions in biofilms, 
in which they are protected against disinfectants. Workers in cooling towers, 
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persons working in air conditioned rooms and people taking a shower are 
exposed to Legionella by inhalation when the systems are not well designed, 
constructed, or maintained. 

Neisseria gonorrhoeae Biofilms 

Neisseria gonorrhoeae is an exclusive human pathogen. Recent studies have 
demonstrated that it utilizes two distinct mechanisms for entry into human 
urethral and cervical epithelial cells involving different bacterial surface ligands 
and host receptors. In addition it has been demonstrated that the gonococcus 
can form biofilms on glass surfaces and over human cells. There is evidence 
for formation of gonococcal biofilms on human cervical epithelial cells during 
natural disease and that outer membrane blebbing by the gonococcus is crucial 
in biofilm formation over human cervical epithelial cells. 

Molecular Genetics 

Technological progress in microscopy, molecular genetics and genome 
analysis has significantly advanced our understanding of the structural 
and molecular aspects of biofilms, especially of extensively studied model 
organisms such as Pseudomonas aeruginosa. Biofilm development can be 
divided into several key steps including attachment, microcolony formation, 
biofilm maturation and dispersion; and in each step bacteria may recruit 
different components and molecules including flagellae, type IV pili, DNA and 
exopolysaccharides. The rapid progress in biofilm research has also unveiled 
several genetic regulation mechanisms implicated in biofilm regulation 
such as quorum sensing and the novel secondary messenger cyclic-di-GMP. 
Understanding the molecular mechanisms of biofilm formation has facilitated 
the exploration of novel strategies to control bacterial biofilms 

Autotroph 

An autotroph, (self-feeding) or producer, is an organism that produces 
complex organic compounds (such as carbohydrates, fats, and proteins) from 
simple substances present in its surroundings, generally using energy from light 
(by photosynthesis) or inorganic chemical reactions (chemosynthesis). They are 
the producers in a food chain, such as plants on land or algae in water. They 
are able to make their own food, and do not need a living energy or carbon 
source. Autotrophs can reduce carbon dioxide (add hydrogen to it) to make 
organic compounds. The reduction of carbon dioxide, a low-energy compound, 
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Fig. 8: Flowchart to determine if a species is autotroph, heterotroph, or a subtype 

creates a store of chemical energy. Most autotrophs use water as the reducing 
agent, but some can use other hydrogen compounds such as hydrogen sulfide. 
Phototrophs, a type of autotroph, convert physical energy from sun light (in 
case of green plants) into chemical energy in the form of reduced carbon. 

Autotrophs can be phototrophs, lithotrophs, or chemotrophs. Phototrophs 
use light as an energy source, while lithotrophs make use of inorganic compounds, 
such as hydrogen sulfide, elemental sulfur, ammonium and ferrous iron, as 
reducing agents for biosynthesis and chemical energy storage. Chemotrophic 
species simply utilize electron donors as a source of energy, whether it be from 
organic or inorganic sources, however in the case of autotrophs, these electron 
donors come from inorganic chemical sources. Phototrophs and lithotrophs 
use a portion of the ATP produced during photosynthesis or the oxidation of 
inorganic compounds to reduce NADP + to NADPH in order to form organic 
compounds. 

HETEROTROPH 

A heterotroph (heteros = "another" "different" and trophe = "nutrition") is 
an organism that cannot fix carbonand uses organic carbon for growth. This 
contrasts with autotrophs, such as plants and algae, which can use energy from 
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Fig. 9: Showing the Overview of cycle between autotrophs and heterotrophs 

sunlight (photoautotrophs) or inorganic compounds (lithoautotrophs) to 
produce organic compounds such as carbohydrates, fats, and proteins from 
inorganic carbon dioxide. These reduced carbon compounds can be used as an 
energy source by the autotroph and provide the energy in food consumed by 
heterotrophs. 

Ecology 

Most heterotrophs are chemoorganoheterotrophs (or simply organotrophs) 
and utilize organic compounds both as a carbon source and an energy source. 
The term "heterotroph" very often refers to chemoorganoheterotrophs. 
Heterotrophs function as consumers in food chains: they obtain organic carbon 
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by eating other heterotrophs or autotrophs. They break down complex organic 
compounds (e.g., carbohydrates, fats, and proteins) produced by autotrophs 
into simpler compounds (e.g., carbohydrates into glucose, fats into fatty acids 
and glycerol, and proteins into amino acids). They release energy by oxidizing 
carbon and hydrogen atoms present in carbohydrates, lipids, and proteins to 
carbon dioxide and water, respectively. 

Most opisthokonts and prokaryotes are heterotrophic; in particular, 
all animals and fungi are heterotrophs. Some animals, such as corals, form 
symbiotic relationships with autotrophs and obtain organic carbon in this 
way. Furthermore, some parasitic plants have also turned fully or partially 
heterotrophic, while carnivorous plants consume animals to augment their 
nitrogen supply while remaining autotrophic. 
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Fig. 10: Stirred bioreactor operated as a 
chemostat, with a continuous inflow (the 
feed) and outflow (the effluent). The rate 
of medium flow is controlled to keep the 
culture volume constant. 



A chemostat (from Chemical 
environment is static) is a bioreactor to 
which fresh medium is continuously added, 
while culture liquid is continuously removed 
to keep the culture volume constant. By 
changing the rate with which medium is 
added to the bioreactor the growth rate of 
themicroorganism can be easily controlled. 

Operation 

Steady State 

One of the most important features of chemostats is that micro-organisms can 
be grown in a physiological steady state. In steady state, growth occurs at a constant 
rate and all culture parameters remain constant (culture volume, dissolved oxygen 
concentration, nutrient and product concentrations, pH, cell density, etc.). In 
addition environmental conditions can be controlled by the experimenter. Micro- 
organisms grown in chemostats naturally strive to steady state: if a low amount of 
cells are present in the bioreactor, the cells can grow at growth rates higher than 
the dilution rate, as growth isn't limited by the addition of the limiting nutrient. 
The limiting nutrient is a nutrient essential for growth, present in the media at 
a limiting concentration (all other nutrients are usually supplied in surplus). 
However, if the cell concentration becomes too high, the amount of cells that are 
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removed from the reactor cannot be replenished by growth as the addition of the 
limiting nutrient is insufficient. This results in a steady state. 
Dilution Rate 

At steady state the specific growth rate (ja) of the micro-organism is equal 
to the dilution rate (D). The dilution rate is defined as the rate of flow of medium 
over the volume of culture in the bioreactor. 

Medium flow rate F 
Culture volume V 
Maximal Growth Rate 

Each microorganism growing on a particular substrate has a maximum 
specific growth rate (|i max ) (the rate of growth observed if none of the nutrients 
are limiting). If a dilution rate is chosen that is higher than |i max , the culture will 
not be able to sustain itself in the bioreactor, and will wash out. Even though 
maximum rates can be obtained, the reactors may become very large. This is 
especially true in E. coli fatty acid production in a glucose medium. 

Applications 

Research 

Chemostats in research are used for investigations in cell biology, as a 
source for large volumes of uniform cells or protein. The chemostat is often 
used to gather steady state data about an organism in order to generate a 
mathematical model relating to its metabolic processes. Chemostats are also 
used as microcosms in ecology and evolutionary biology. In the one case, 
mutation/selection is a nuisance, in the other case, it is the desired process 
under study. Chemostats can also be used to enrich for specific types of bacterial 
mutants in culture such as auxotrophs or those that are resistant to antibiotics or 
bacteriophages for further scientific study. Competition for single and multiple 
resources, the evolution of resource acquisition and utilization pathways, cross- 
feeding/symbiosis, antagonism, predation, and competition among predators 
have all been studied in ecology and evolutionary biology using chemostats. 
Industry 

Chemostats are frequently used in the industrial manufacturing of ethanol. 
In this case, several chemostats are used in series, each maintained at decreasing 
sugar concentrations. 
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Concerns 



Foaming results in overflow with the volume of liquid not exactly 

constant. 
Some very fragile cells are ruptured during agitation and aeration. 
Cells may grow on the walls or adhere to other surfaces, which is easily 

overcome by treating the glass walls of the vessel with a silane to 

render them hydrophobic. 
Mixing may not truly be uniform, upsetting the "static" property of the 

chemostat. 
Dripping the media into the chamber actually results in small pulses of 

nutrients and thus oscillations in concentrations, again upsetting the 

"static" property of the chemostat. 
Bacteria travel upstream quite easily. They will reach the reservoir of 

sterile medium quickly unless the liquid path is interrupted by an air 

break in which the medium falls in drops through air. 

Continuous efforts to remedy each defect lead to variations on the basic 
chemostat quite regularly. Examples in the literature are numerous. 

Antifoaming agents are used to suppress foaming. 

Agitation and aeration can be done gently. 

Many approaches have been taken to reduce wall growth. 

Various applications use paddles, bubbling, or other mechanisms for 

mixing. 
Dripping can be made less drastic with smaller droplets and larger vessel 

volumes. 
Many improvements target the threat of contamination. 
Variations 

Fermentation setups closely related to the chemostats are the turbidostat, 
the auxostat and the retentostat. In retentostats culture liquid is also removed 
from the bioreactor, but a filter retains the biomass. In this case, the biomass 
concentration increases until the nutrient requirement for biomass maintenance 
has become equal to the amount of limiting nutrient that can be consumed. 
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TURBIDOSTAT 

A turbidostat is a continuous culture device, similar to a chemostat or an 
auxostat, which has feedback between the turbidity of the culture vessel and the 
dilution rate. The theoretical relationship between growth in a chemostat and 
growth in a turbidostat is somewhat complex, in part because it is similar. A 
chemostat technically has a fixed volume and flow rate - thus a fixed dilution rate. 
When the cells are uniform and at equilibrium, operation of a chemostat and 
turbidostat should be identical. It is only when classical chemostat assumptions 
are violated (for instance, out of equilibrium; or the cells are mutating) that a 
turbidostat is functionally different. One case may be while cells are growing at 
their maximum growth rate, in which case it is difficult to set a chemostat to the 
appropriate constant dilution rate. 

While most turbidostats use a spectrophotometer/turbidometer to 
measure the optical density for control purposes, there exist other options, such 
as dielectric permittivity. 

A turbidostat uses a continuous culturing method where the turbidity of 
the culture is held constant by manipulating the rate at which medium is fed. If 
the turbidity tends to increase, the feed rate is increased to dilute the turbidity 
back to its set-point. When the turbidity tends to fall, the feed rate is lowered so 
that growth can restore the turbidity to its set-point. 

The problem of growth or other materials fouling the optical surfaces of 
whatever method is used to measure turbidity has not been solved. While a 
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Fig. 11: Sketch of Turbidostat 
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turbidostat may operate well for a brief time, the control signal for turbidity 
soon becomes unreliable. 

AUXOSTAT 

A device that uses the rate of feeding to control a state variable in continuous 
culture is termed an auxostat. The organisms establish their own dilution 
rate. While the well-known chemostat is stable and simple for investigating 
continuous cultivation at low to moderate dilution rates, an auxostat tends to 
be much more stable at high dilution rates. Population selection pressures in 
an auxostat lead to cultures that grow rapidly. Practical applications include 
high-rate propagation, destruction of wastes with control at a concentration 
for maximum rate, open culturing because potential contaminating organisms 
cannot adapt before washing out, and operation of processes that benefit from 
careful balance of the ratios of nutrient concentrations. 

This name was coined by Martin and Hempfling although earlier 
investigators had proposed the terms nutristat, nustat, and controlled- 
concentration-coupled-continuous-cultivation (C5). 
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Fig. 12: Showing the Schematic Diagram of the Auxostat Design and Operation 
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Some Calculations of Auxostat Operation 



Feedback control of nutrient or product concentrations in bioprocesses 
has been progressing in fits and starts for 40 years (Bryson, 1952), but has 
suffered from the lack of reliable equipment. Continuous bioprocessing permits 
the observation of microorganisms at steady state. Although chemostats are in 
common use, the advantages of auxostats have been recognized for sometime 
(Agrawal and Lim, 1984; Fraleigh et al, 1989). Auxostats can be operated in 
difficult or unstable conditions and can achieve steady state more quickly than 
can the open-loop chemostat. Operation near the maximum growth rate can 
be extremely unstable with a chemostat. The following equations are based on 
a bioprocess having Monod kinetics with mass balances for the liquid phase of 
the bioreactor: 



x- BY 



D^MY 



dP 
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where iJL = specific growth rate coeff. 

it - maximum growth rate coeff. 

S = concentration of limiting nutrient 
K s = half-saturation coeff. 

nsc-ns= y + mx 

X = biomass concentration (g l" 1 ) 

P = product concentration (g l" 1 ) 

D = dilution rate (h" 1 ) 

Yis usually Y x / S , biomass yield on nutrient (g g" 1 ) 

Other important yield coefficients are: 

Y / x = product yield on biomass (g g" 1 ) 

Y y s = product yield on substrate (g g" 1 ) 
V = volumetric productivity (g l" 1 ) 
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Operation of an Auxostat 



At steady state the dilution rate is equivalent to [i. In a chemostat the dilution 
rate is the selected feed rate divided by the working volume, and the biological 
parameters must adjust to it, Le., open-loop operation. At high growth rates, 
the values of the biological state variables (X, S, Vp) are very sensitive to small 
changes in dilution rate. In the terminology of control engineers, this is a high 
gain region of operation in which a small change in dilution rate relates to large 
changes in concentrations of cells and nutrient with the well-known, inherent 
problem of "washout". If the dilution rate rises above the maximum growth rate 
of the microorganism, the effective biomass concentration becomes zero. 

Fluctuations in dilution rate may result from drift in pumping rates and 
changes in the bioreactor working volume (D = flow rate/volume). In bench 
scale bioprocessing, peristaltic pumps are susceptible to short and long term 
changes in flow rates. The working volume is normally assumed to be constant 
and is rarely measured on-line, but it can fluctuate significantly with changes 
in aeration, surface tension, and biomass concentration. Therefore, the desired 
steady state in continuous culture is difficult to achieve with a chemostat as the 
dilution rate approaches the maximum specific growth rate, and this is often the 
region of highest productivity. 

Auxostats eliminate the operational problems found at high growth rates. 
Stability arises in the bioreactor through manipulation of the pumping rate of 
fresh medium using feedback control strategies. The dilution rate becomes 
a dependent variable that adjusts to the rates of supply and demand for the 
controlled variable. Random short-term disturbances in the pumping rate 
and bioreactor volume are handled by the control algorithm and long-term 
changes can be accounted for by volume measurements and intermittent 
pump calibration. This permits stable operation in the 'high gain areas near the 
maximum growth rate, but only when robust and accurate measurement and 
control for the state variables can be devised. Fig. 12 is a sketch of an auxostat. 

Uses and Advantages of Auxostat 

Auxostats have an inherent advantage because they can reach steady state 
more rapidly than can chemostats. Chemostats have a very poor response time 
as the dilution rate approaches the maximum growth rate, where it increases 
exponentially. Faster response times are a basic benefit of feedback control 
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(Seborg et ah, 1989). Disturbances and changes in setpoints are very common 
in bench scale bioprocessing (planned and unplanned), and feedback control 
permits quicker attainment of a new steady state. This allows auxostats to 
move quickly through different operational regimes. The choice of the feedback 
variable for an auxostat is quite broad, however it must relate to microbial growth 
(Fredrickson etaL, 1970). Some of the earliest work employed turbidostats where 
the biomass concentration, as reflected by turbidity, was controlled through the 
manipulation of dilution rate (Bryson, 1952). This technique is normally based 
on measuring optical density, however it is unreliable because of fouling of the 
sensing element (Watson, 1969). It is still used as research tool for non-fouling 
types of cultures (Aarinio etal, 1991). 

Types of Auxostats 

Dissolved Oxygen Auxostat 

Product and nutrient based auxostats use direct measurements of chemical 
species that directly reflect metabolism. The types are limited because of the 
lack of suitable sensors. Control based on the dissolved oxygen concentration 
(DO) was first proposed by Ohashi (1958). Operation is based on manipulating 
the liquid feed rate to maintain a certain DO level and not the manipulation of 
the air flow rate. Even though oxygen is a nutrient, it is normally not limiting in 
this type of operation and the absolute value is relatively unimportant. The DO- 
auxostat is similar in nature to a product-based auxostat because the dilution 
rate does not have a direct effect on DO. The flow of fresh air is not manipulated 
in simple operations but contributes along with the agitation rate to an overall 
mass transfer coefficient (Kja), that is assumed to be constant. The governing 
equation for this type of auxostat is shown in the equation, that is based on an 
2 mass balance. 

^-= RX + K|a[C*-C) 

where C = concentration of oxygen 

C*= concentration at saturation 
K|a = lumped oxygen transfer concentration 
R = uptake rate of oxygen by organisms 
X = concentration of organisms 




The term in the equation that has the coefficient times the concentration 
difference is generally termed the oxygen transfer rate (OTR) and the term with 
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the minus sign is the oxygen uptake rate (OUR), By manipulating the feed rate of 
a limiting nutrient such as sugar, the OUR can be changed thereby effecting the 
DO. Oxygen consumption must be growth-related for effective control. 

The first control of dissolved oxygen was reported with a yeast bioprocess 
using on/off control of the feed rate (Hospodka, 1966). The time dependency of 
the DO concentration and dilution rate were not reported. Yamada, etal (1979) 
used this technique for controlling the conversion of sorbitol to sorbose by 
Acetobacter suboxydans. They also implemented control based on the oxygen 
content of the bioreactor off-gas and obtained a more stable operation than 
with dissolved oxygen concentration. To improve the control, Konstantinov, 
et al (1990) used a more advanced technique where unmeasured, small, rapid 
disturbances were rejected through manipulation of the agitation rate, which 
adjusted the OTR. These disturbances were attributed to fluctuations in air flow 
rate and pressure. Feed rate was changed when the agitation rate moved outside 
a predetermined window of operation to adjust for OUR changes. 

Dual Setpoint Auxostat 

Introduction 

Continuous bioprocessing permits the observation of microorganisms 
at steady state. Traditionally chemostats are used, however, the advantages of 
auxostats in which a concentration in the bioreactor is held constant by varying 
the nutrient feed rate have been recognized for sometime (Agrawal and Lim, 
1984, Fraleigh, et aL, 1989). An auxostat is generally more stable than the 
open-loop chemostat at dilution rates near the maximum growth rate of the 
microorganism. Instabilities arise because of fluctuations in pumping rate and 
reactor volume, but in an auxostat the dilution rate becomes a dependent variable 
and is a function of the setpoint of the controlled variable. Chemostat operation 
with toxic nutrients or chemicals is impossible at inhibiting concentrations, but 
auxostat stability is achieved through manipulation of the pumping rate of fresh 
medium using feedback control strategies to maintain a desired concentration. 
pH Auxostat 

The pH-auxostat couples the addition of fresh medium to pH control. As 
the pH drifts from a given setpoint, fresh medium is added to bring the pH back 
to the setpoint. The rate of medium addition is determined by the buffering 
capacity and the feed concentration of the limiting nutrient and not directly 
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by the setpoint (pH) as in a traditional auxostat. The pH-auxostat is robust but 
controls nutrient concentration indirectly. The pH change is often an excellent 
indication of growth and meets the requirements as a growth-dependent 
parameter as defined by Fredrickson etal (1970). However, the exact cause of pH 
change varies among organisms. It represents the summation of the production 
of different ionic species and ion release during substrate uptake. Therefore the 
pH can move either up or down as a function of growth. The most common 
situation is pH depression because of organic acid production and ammonium 
uptake. However for microorganisms growing on protein or amino acid-rich 
media, the pH will rise with growth because of the release of excess ammonia. 
The pH-auxostat was simultaneously reported by different research groups and 
was probably first operated in the 1960s (Bungay, 1972; Watson 1969) and a 
detailed development was presented by Martin and Hempfling (1976). 
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Fig. 13: Showing the Lay out Diagram of a pH Auxostat 
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The pH-auxostat is the most widely applied type of auxostat. Oltmann 
et al. (1978) evaluated it for efficient, high-productivity biomass production. 
Their improvement was separating the buffer from the nutrient medium to 
overcome inhibition of high buffer concentrations. Fraleigh (1988) used a pH- 
auxostat to study aerobic production of ethanol by yeast operating near the 
maximum specific growth rate where specific ethanol productivity was highest. 
It also permitted an investigation of the coupling of energy and acid production 
that was not possible in a chemostat. 

Using two reservoirs at different pH expedites changing the buffer capacity 
simply by changing the proportions. Selecting new pumping rates changes the 
buffering capacity. 
Controllers of Auxostats 

The direct and independent control of two or more nutrient and/or product 
concentrations by using separate media streams should open new avenues 
of research. The ratio of the limiting nutrient to other media components is 
critical in the production of many bioproducts such as storage compounds, 
exo-polysaccharides, and enzymes (Dawes and Senior, 1973, Kole, {et al.}, 
1988). Multi-substrate limited growth, where more than one nutrient effects the 
growth rate is also an important area that has been neglected (Egli and Quayle, 
1986, Rutgers, {et al.}, 1990). These types of interactions can be explored with a 
chemostat by manipulating the medium concentrations, but this is cumbersome 
and indirect. True multivariable control has two major strengths: 

1. direct control of the nutrient levels in the fermentor, 

2. easy and quick manipulation of nutrient ratios. 

The ability to set the nutrient concentrations directly is a significant 
improvement over the chemostat. In a chemostat, a nutrient ratio is a dependent 
variable, subject to many influences that make reproducibility extremely 
difficult. Feedback control also improves response times, allowing quicker 
identification of ideal operating conditions. In addition, the response dynamics 
of setpoint changes should give insight into the interaction between substrate 
uptake and metabolism. 

Multivariable control in continuous culture is well-known but is 
generally limited to variables such as pH, dissolved oxygen concentration, and 
temperature. In these cases, the parameter is only indirectly related to growth, 
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the medium feed rate is not the manipulated variable, and the control loops are 
independent. Multivariable control of nutrient concentrations directly related to 
growth and metabolism has been implemented in fed-batch cultures. Kole etal. 
(1988) improved protease production by B. subtilis using multivariable control. 
Ammonium was controlled directly and glucose indirectly by adding glucose- 
containing media as a function of the dissolved oxygen level. Suzuki etal. (1986) 
investigated production of poly-beta-hydroxybutyrate in fed-batch cultures of 
a methylotrophic Protomonas species. Methanol levels were controlled directly 
through the entire fermentation and ammonium indirectly during the growth 
phase by using ammonia solutions for pH control. They determined that the 
initial growth phase was shortened by this multivariable control technique. 

Ethanol Auxostat 

An ethanol-auxostat was operated in an acetic acid bioprocess using an 
organic vapor sensor (Muhlemann and Bungay, 1993). Ethanol concentration 
in a bioreactor was held constant by using the signal from a commercial gas- 
sensor for organic vapors. Gas-sensors for alarm systems are very inexpensive 
and are suitable for organic nutrients that are volatile. Unlike an electrode, a 
gas-sensor is electrically insulated from the liquid environment and is unlikely 
to become fouled. Substrate concentration in the gas phase above the broth 
should be in dynamic equilibrium with the liquid phase. 

This gas-sensor consists of a sensing element and a heater. The small sensing 
element is Sn0 2 sintered metal. When heated to about 400 C with no oxygen 
present, free electrons flow easily through the grain boundary of the tin dioxide 
Sn0 2 particles. Oxygen adsorbed on the particle surface forms a potential 
barrier in the grain boundaries that traps free electrons. This restricts their flow, 
causing the electrical resistance to increase. When the sensor element adsorbs 
reducing gases, their oxidation lowers the potential barrier, allowing electrons 
to flow more easily, thereby reducing the electric resistance. 

The sensor mounted in the fermenter head space can be seen in the next 
figure. The unit in a rubber stopper fits into the fermenter lid. Air from the head 
space flows through the flame arrester to the sensing element. To get renewal 
around the sensing element, gas is allowed to escape through the back of the 
sensor. A stainless steel grid shaped as a cone prevents splashing from reaching 
the sensor. 
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Fig. 14: Showing the Diagram of the Ethanol Auxostat 

For aseptic operation, the electrode should be mounted in the exit gas line. 
The additional lag time would be small. Calibration of the whole system links the 
ethanol concentration in the broth with the sensor output voltage. Air stripping 
of ethanol from the culture fluid was measured and could be reconciled in mass 
balances. 

An acetic acid bioprocess was the model for evaluating the ethanol auxostat. 
Ethanol was the growth-limiting substrate. Soil was the inoculum. The auxostat 
principle was demonstrated by runs of several hundred hours. All bioreactors 
were open, and only the feed was sterile. Use "BACK" to go with the flow. 

Analysis of Buffering Capacity 

Buffering capacity is defined as the equivalents of titrant required to change 
the medium pH to the reactor pH. The additional governing equation is the 
mass balance on the H + ion concentration in the bioreactor. 

By using fresh medium for pH control, the limiting nutrient level in the 
reactor can be manipulated through adjustments in buffering capacity and/or 
limiting nutrient concentration in the feed. Assuming steady state and that the 
difference between feed and bioreactor H + is very small (pH 7 means a very tiny 
amount of H + ), simplification gives the following: 
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An extensive study of buffering capacity by Rice and Hempfling (1985) 
showed that the specific growth rate was roughly ten times greater with low 
buffering capacity than with a well-buffered system. The pH-auxostat is the most 
widely applied type of auxostat. Oltmann, etaL, (1978) evaluated it for efficient 
high-rate biomass production. Their improvement was separating the buffer 
from the nutrient medium to overcome inhibition of high buffer concentrations. 
Fraleigh, (1988) used a pH-auxostat to study aerobic production of ethanol by 
yeast operating near the maximum specific growth rate where specific ethanol 
productivity was highest. It also permitted an investigation of the coupling of 
energy and acid production, which was not possible in a chemostat. 

System Crashes 

Minkevich, etal., (1988) improved the pH-auxostat to overcome buffering 
capacity mismatch. If the buffering capacity is set too high, there will be 
insufficient substrate to maintain the biomass concentration determined by the 
buffering capacity. In this situation the pH will not drop below the set-point 
(assuming that acids are produced) before depleting the reactor of the limiting 
nutrient. Therefore the media pumps will not turn on because without nutrient 
to metabolize, the pH will not change. To overcome this problem, a fixed feed 
rate component should be added to the overall feed rate. 
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